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IDENTIFICATION AND MEASUREMENT OF COMBUSTION NOISE FROM A TURBOFAN 
ENGINE USING CORRELATION AND COHERENCE TECHNIQUES 


Abstract 

by 

ALLEN MARTIN KARCHMER 

A 


Fluctuating pressure measurements within the combustor and tail- 
pipe of a turbofan engine are made simultaneously with far-field 

* 

acoustic measurements. The pressure measurements within the engine 
are accomplished with cooled '’semi-infinite" waveguide probes utiliz- 
ing conventional condenser microphones as the transducers. The mea- 
surements are taken over a broad range of engine operating conditions 
and for sixteen far-field microphone positions between 10° and 160° 
relative to the engine inlet axis. Correlation and coherence tech- 
niques are used to determine the relative phase and amplitude re- 
lationships between the internal pressures and far-field acoustic 
pressures. The results indicate that the combustor is a low fre- 
quency source region for acoustic propagation through the tailpipe 
and out to the far-field. Specifically, it is found that the rela- 
tion between source pressure and the resulting sound pressure in- 
volves a 180° phase shift. This latter result is obtained by 
Fourier transforming the cross-correlation function between the 
source pressure and acoustic pressure after removing the propagation 
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delay time. Further, it is found that the transfer function between 
the source pressure and acoustic pressure has a magnitude approxi- 
mately proportional to frequency squared. These results are shown 
to be consistent with a model using a modified source term in Light- 
hill * s turbulence stress tensor, wherein the fluctuating Reynolds 
stresses are replaced with the pressure fluctuations due to fluc- 
tuating entropy. Finally, the measured ordinary coherence functions 
between fluctuating combustor pressure and far field acoustic pres- 
sures are used to obtain the far-field combustor associated noise in 
terms of spectra, directivity and acoustic power. The results so 
obtained are compared with various existing empirical prediction pro- 
cedures and are found to be in good agreement with respect to direc- 
tivity patterns, but somewhat low when compared with respect to 


acoustic power. 
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1. INTRODUCTION 


Noise, defined as unwanted or excessive sound, is a relatively 
local phenomenon affecting only the immediate vicinity of its 
source, and does not accumulate in the environment , Nevertheless 
it is now broadly recognized as a form of environmental degrada- 
tion, Noise can be annoying, can interfere with sleep, work, or 
recreation and, under some circumstances, may cause physical and 
psychological damage. 

The most prevalent noise sources in society today stem from 
transportation operations with aircraft noise playing a particu- 
larly significant role, especially since the introduction of jet 
powered aircraft into commercial service in the late 1950 f s* Be- 
cause of increasing public concern with the environment and the 
phenomenal growth of commercial aviation in this country (predicted 
to reach 430 million flight operations per year by 19? S (ref, 1)) 
aircraft noise abatement has become a major objective of both in- 
dustry and government. 

Efforts to reduce airport community noise have essentially 
taken two paths* The first is through a multitude of aircraft and 
airport operational procedures such as the multi-segment approach, 
power cutback and turnouts after takeoff, runway restrictions, the 
establishment of airport buffer zones, land use control, and even 
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curfews. Operational restrictions such as these, however, can only 
be partially effective or are limited in their implementation be- 
cause of local conditions or safety considerations. The second ap- 
proach, and the one where the preponderance of research effort is 
concentrated is aimed at the reduction or alteration of the primary 
aircraft noise source itself, the propulsion system. 

A jet powered propulsion system is a complex acoustic source. 
It is composed of a number of individual but interrelated noise 
sources. While these sources are not strictly independent of one 
another, we can generally associate with each of the noise sources 
within a turbofan engine some identifiable physical generating 
mechanism and acoustic characteristics which, in turn, are usually 
related to specific components of the engine. It is entirely ap- 
propriate, therefore, when analysing the constituents of aircraft 
propulsion system noise to analyse them (or synthesize, when pre- 
dicting) on an engine component by component basis. 

Figure 1, for example, is an illustration of a typical turbo- 
fan engine, schematically depicting the primary component noise 
sources. The noise which is emitted from the inlet and from the 
fan discharge duct is generated principally by the fan and, to a 
lesser extent, by the compressor. The fan is responsible for the 
high frequency whine heard principally on approach and primarily 
forward of the aircraft. The noise emitted by the primary jet ex- 
haust at the rear of the engine is generated by the violent turbu- 
lent mixing between the high velocity exhaust jet and the rela- 



tively quiescent ambient environment. This turbulent mixing noise 
is responsible for the low frequency roar which is heard principally 
on takeoff and primarily in the rear of the aircraft. Less impor- 
tant, but not insignificant, is the noise generated by the second- 
ary, or fan, jet exhaust as well as turbine noises which are emitted 
at the rear of the engine. Finally, there are noises associated 
with the combustion of the air-fuel mixture, as well as noises due 
to the interaction of the engine gases as they flow through, over 
and around various internal surfaces and passages. 

These latter sources, because they originate in the gas gen- 
erator section, or core, of the turbofan, have collectively been 
given the name core noise. For virtually all of the aircraft in 
commercial service today, fan and jet noise are thought to be, by 
significant margins, the dominant contributors to overall engine 
noise, with core noise sources relatively unimportant. 

In recent years, however, there have been substantially in- 
creased research efforts to determine the causes, generating mech- 
anisms, and characteristics of core associated noise in turbofan 
engines. There are a number of reasons for this and these will be 
briefly discussed. 

The general conclusion that the fan and jet are the primary 
sources of turbofan engine noise has, for the most part, been 
reached as a result of extensive model nozzle, fan and full scale 
engine static tests; that is, tests under nonflight conditions. It 
has long been recognized, however, that jet mixing noise is signif- 
icantly reduced by forward motion as a result of the reduction of 


the relative velocity between the jet and the ambient atmosphere. 

One of the earliest analyses of jet noise (by Lighthill, ref. 2) in- 
dicated that the total sound power generated by a turbulent jet is 
proportional to the eighth power of jet velocity. This has since 
been verified by numerous investigators. Relatively modest reduc- 
tions in jet velocity s therefore, provide enormous reductions in jet 
noise. Reduction of jet velocity by 50 percent, for example, results 
in a 250 fold (24 dB) reduction in sound power. Model jet experi- 
ments in various forward velocity simulation facilities (e.g., large 
free jets, wind tunnels, etc.) have shown similar strong dependen- 
cies by the jet noise on the relative velocity between jet and ambi- 
ent flow, although not as strong as eighth power (refs. 3, 4, 5). 
Forward motion reduces relative velocity and therefore reduces jet 
noise generation. 

Recent experimental studies (refs. 6 and 7, for example) have 
also shown reductions in fan noise generation under simulated for- 
ward velocity conditions, primarily as a result of inflow cleanup 
conditions. That is, the forward velocity conditions eliminate or 
minimize the ingestion of ground vortices or large scale turbulence, 
thus reducing the noise generated as these flow Irregularities pass 
through the fan. 

A significant theoretical effort has also been made in attempts 
to understand such forward flight effects on fan and jet noise (refs. 
3 and 6). For the most part, these analytical approaches are able to 
predict the acoustic characteristics observed on model nozzles and 
fans in simulation facilities. However, a comparison of the theo- 


retically predicted or experimentally observed model simulation re- 
sults with the noise characteristics measured from actual aircraft 
flight tests shows large discrepancies. The flight data generally 
indicate more noise than would be expected from the theoretical pre- 
dictions or model simulation tests. The reasons for these discrep- 
ancies are not fully understood. Some investigators (refs. 8 and 9), 
however, have advanced the argument that internal noise sources 
(i.e., core noise), which are not simulated in model nozzle or fan 
tests, and which are not observable in static engine tests because 
of extremely high jet noise levels, become significant contributors 
to overall engine noise during takeoff or approach when the fan and 
jet noise are reduced because of forward velocity effects. This 
would indicate that core noise, consequently, may present a floor 
to the aircraft noise abatement engineer. If new federal regula- 
tions mandate further engine noise reductions, it is possible that 
the effort will have to be aimed at these core noise sources as well 
as the fan and the jet. 

The second factor stimulating research into core noise is re- 
lated to the possible development of an advanced supersonic cruise 
commercial aircraft. In the early 1970's it became obvious that with 
the then existing technology, a commercial supersonic aircraft would 
generate unacceptable levels of jet noise because of the extremely 
high jet velocities. At least partially for this reason, further 
development work on an American SST was halted. Continued research, 
however, has recently indicated that inverting the conventional tur- 
bofan jet exhaust velocity profile so that the outside (fan) jet has 
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a greater velocity than the inside (primary, or core) jet results in 
a significant reduction in jet noise without correspondingly signif- 
icant thrust losses (refs. 10 and 11). 

One proposed method of accomplishing this inversion simply re- 
quires direct supplementary combustion in the fan duct. In the con- 
ventional turbofan, it is generally thought that the presence of the 
turbine between combustor and exhaust significantly attenuates com- 
bustion associated noise (between 5 and 20 dB, depending on fre- 
quency, see ref. 12). It is possible, therefore, that for the duct- 
burning turbofan, where there is no intervening turbomachinery down- 
stream of the supplementary combustor, combustion noise may be a 
major contributor to core noise. 

It should be noted that although the terms "combustion" noise 
and "core" noise are often used interchangeably in the literature * 
they refer to distinctly different phenomena. In the present work, 
the tern "core" noise refers collectively to all noise or noise 
sources associated with the gas generator of a turbofan. As noted 
previously this can include compressor noise, turbine noise, noises 
associated with the gases flowing through, over and around the vari- 
ous internal surfaces and passages, as well as combustion noise. 

The term "combustion noise" refers only to the noise generated as a 
result of the combustion process (direct or indirect), and which 
would be absent if combustion did not take place. This is not a 
trivial distinction and, in fact, provides a primary motivation for 


the present work. 


A typical procedure in identifying core noise from an engine is 
to make a series of far-field acoustic measurements at progressively 
lower engine power settings. Fan, compressor, and turbine noise can 
generally be Identified in the data because of the characteristic 
tone-like spectra associated with their high rotational speeds. (Al- 
ternatively, the fan noise can be muffled, as in ref. 13, and the jet 
noise can be reduced by replacing the nozzle with a section of 
straight pipe, as in ref, 14.) The data are then examined for be- 
havior not characteristic of jet noise (e.g., sound power not pro- 
portional to the eighth power of jet velocity; improper Strouhal 
scaling of spectra, etc.). The observed differences are then at- 
tributed to internal, or core, noise sources. Such a technique, at 
best, is imprecise from a quantitative viewpoint, and is unable to 
qualitatively distinguish one internal source from another. 

Several investigators (refs. 15 and 16) have attempted to 
identify specific regions within the core as contributors to far- 
field noise by combining internal pressure measurements with far- 
field acoustic measurements and cross-correlating them. While this 
approach is generally more useful than acoustic measurements taken 
only outside the engine, cross-correlation measurements by themselves 
shed very little light on cause-effect relationships. For example, 
it is certainly possible, as will be shown later, to have a non-zero 
correlation between the pressure (or any other fluid property) at 
some point within the engine and the acoustic pressure outside the 
engine and yet have the source of the acoustic information at some 
unidentified other point within the engine. That is, correlation 
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alone does not demonstrate cause-effect, although examination of time 
delay information can provide useful information, as was shown in 
reference 15. 

There have also been numerous attempts using a variety of mea- 
surement techniques to specifically identify and characterize the 
noise from gas turbine combustors in non-engine installations. 

That is, such investigations are conducted on combustors installed 
in specially designed apparatus which do not contain the addi- 
tional machinery and hardware (compressor, turbine, nozzle, etc.) 
associated with a turbofan propulsion system. The techniques range 
from correlation analysis (ref. 17) and coherence analysis (ref. 18) 
to simple far field acoustic measurements (ref. 19). In experiments 
such as these the noise source, of course, is known and the effort 
involves one of quantifying the combustion noise contribution. It 
is not clear, however, how the results of these experiments can be 
applied or related to an actual engine where the presence of the 
turbine and nozzle can serve to partially attenuate or alter some 
of the noise generated by the combustor, as well as generate noise 
themselves. 

The present work will describe an experiment conducted on an 
operating turbofan engine. The objective of the investigation was 
to determine specifically whether or not the combustor could be iden- 
tified as a contributor to far-field engine noise and to quantify 
this contribution. The engine was not modified in any way to reduce 
the noise generated by the fan or jet. Correlation, coherence, and 
cross-spectral measurements were made between fluctuating pressures 
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at various locations within the gas-generating section of the engine 
and acoustic pressure in the far-field of the engine. 

The results of the measurement program are basically divided 
into two parts: a diagnostic phase which focuses on a rigorous in- 

terpretation of certain correlation function characteristics to de- 
termine the relationship between internal pressures and far-field 
acoustic pressure and thus identify acoustic source regions and pos- 
sible physical generating mechanisms; and a computational phase 
which presents the quantitative contribution of the combustor as- 
sociated far-field noise to overall engine noise in terms of spectra, 
directivity, and level, including acoustic power. 


t I 1 I l T 1 


2, PHYSICAL CONSIDERATIONS 
2.1 Derivation of Ligh thill’s Equation 

As indicated in the INTRODUCTION , the existence of a correla- 
tion between two physical variables in no way implies a cause-effect 
relationship. However, there is additional information on the rela- 
tive phase and amplitude between the two variables being correlated. 
This information can be extracted through the frequency domain rep- 
resentation of the correlation function by the appropriate Fourier 
transformation of the correlation function. It is this phase and 
amplitude information which can provide clues as to the nature of 
the relationship between the two variables. 

For example, one expects that if the pressure in an acoustic 
source region were correlated with the pressure at an acoustic field 
(i.e,, observation) point, then the phase and amplitude relationship 
derived from that correlation should be consistent with some physical 
generating, or source mechanism. If, on the other hand, the two 
pressures being correlated were both simply in the acoustic field of 
a third point which is a common source for those two, then the phase 
and amplitude relationship should be consistent with pure time de- 
lay. It is this notion upon which the diagnostic phase of the pres- 
ent work is based. It is appropriate at this point, therefore, to 
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briefly review the derivation and interpretation of the governing 
equations for aerodynamic noise generation. 

In 1952 and 1954, in two now famous papers (refs. 2 and 20), 

M. J. Lighthill proposed his fundamental theory of noise production 
by aerodynamic sources alone. Physically, he showed that sound is 
generated by an unsteady fluid flow exactly as in a uniform medium 
at rest which is acted upon by certain externally applied fluctuating 
stresses. His derivation is as follows. 

The starting point is the continuity equation 


9p 

3t 


3pu. 


oXj 


- 0 


and the momentum equation (neglecting body forces) : 

Du. 9e' . 

i = ijL 

P Dt 3Xj 


where 


( 1 ) 


( 2 ) 


_D_ _ _9_ , 9 

Dt _ 9t U j 9x_. 

and e ! . is the stress tensor, 
ij 

The equations are written for an infinite, uniform quiescent 
medium, which contains a region of turbulent flow. Here p = den- 
sity, and u^ is the component of velocity in the x^ direction 
(i - 1,2,3). 

Expanding the material derivative in (2), adding u^(3pu./9Xj) 
to both sides, and substituting continuity into the result gives: 
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(3) 


Differentiating continuity with respect to time, and differentiating 
(3) with respect to x^ gives, respectively: 


a 2 n 9 PU i 

JLfi. + — = o 

3t 2 St 3x i 
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and 
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3t 9Xj ax^^ ^ 6 ij PU i U j^ 


(5) 


2 2 

Subtracting (5) from (4) and adding t^(8 p/Sx^ 3x^) to the result 
gives : 


c 2 .-^g— - ■ ■■ 2 . ' lou.u. - e;, - c‘p«„| (6) 

8t 2 0 3x i 3x i 3x i 3x j 


: frVj - e ij ■ c2si d 


where is the constant speed of sound and is the Kronecker 

Delta, 


For an isotropic Newtonian fluid* 


e lj = _6 ijP + 2ve ij + 6 ij Xe ii 


(7) 


where 


p is the shear coefficient of viscosity 
X is the bulk coefficient of viscosity 

is the strain rate tensoif = 1/2 [ (3u^/3xj) + (3Uj/3x^)] 
p is the pressure 

For the traditional assumption of a Stokesian fluid, X + 2/3 p = 0 
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and (6) becomes: 


2 9 2 

i_a 2 r 0 = r 

^,.2 0 3x. 9x. ~ 9x. 9x 

9t x x X j 


: [ pu i u 3 - e ij + (p - 


C8) 


where e. . is the viscous part of the stress tensor: 
ij 


e. . 
3-1 


3u^ 9u . 

9x. + 9x. 

3 i 


3 3^ °ij 


Finally, then, equation (8) becomes: 


2 2 
lie. _ r 2 -IIP- 
2 

9t 


0 9x. 9x. 
xi 


9x. 9x. (T ij 5 

-t -i *■* 


where 


T u = pu i “ 3 - •« + 'p - c p)S ij 
Referring the pressure and density to (constant) reference 
values (subscript 0) : 

p’ = P - Pn 


P T “ P - P, 


Then 


l?£l _ c 2 aV-. = ... 3 2 - - (T , J 

at 2 0 9 Xi 9x i 3x ± 9sc u ij' 


(9) 


where T!. is known as Lighthill’s turbulence stress tensor, 
xj 

Equation (9) clearly has the form of a wave equation describing 
the propagation of sound in a stationary medium driven by a distribu- 
tion of quadrupole sources 9^/3x. 3x.(T!.). This analogy between 
the density fluctuations in a real flow and those in an ideal acous- 
tic medium at rest embedded with a distribution of quadrupoles is 
called Lighthill T s acoustic analogy. 
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More importantly, however, is the fact that equation (9) is an 
exact equation, at least within the assumptions of the momentum equa- 
tion, continuity equation, and the assumption of a Newtonian and 
Stokesian medium. All real flow and acoustic effects are accounted 
for: refraction, convection, turbulence scattering of acoustic 

waves, mean flow-acoustic interaction effects, and viscous dissipa- 
tion of the sound waves by the flow. 

Some difficulty arises upon a strict interpretation that the 
right hand side is a pure source term, and the left-hand side ac- 
counts solely for acoustic propagation. The source term, , con- 
tains the fluctuating density, which is the dependent variable 
being solved for on the left-hand side. Since T! . is not, in gen- 
eral, an a priori known quantity except for the most trivial of 
flows, a complete description of the "source term" would imply a 
solution to the general non-linear flow equations. Note that all 
the non-linear terms are on the right hand side of (8) . The precise 
description of the effects of refraction, convection, acousti'c-flow 
interactions, etc. would require that a solution of (8) already be 
in hand. 

Nevertheless, Lighthill's equation has been extensively used 
over the years, particularly for jet noise analysis, where the lack 
of solid surfaces in the region under consideration permits rela- 
tively simple boundary conditions in its solution. Various approach- 
es to incorporate more of the real fluid effects in the wave operator 
portion of the equation and less in the source term have been de- 
veloped (for example Lilley in ref. 21, Phillips in ref. 22), These 
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procedures necessarily deviate from the acoustic analogy approach 
of Lighthill and although more direct and precise from a computa- 
tional viewpoint, are less intuitively satisfying from a physical 
viewpoint . 

2.2 Solution of Lighthill ' b Equation 


Lighthill obtained the solution to his equation, in the absence 
of solid surfaces, as: 


p'(X,t) = p(X,t) - 



rH 

B^CY.t’) 

|X- Y( 

^i^Yj 


d^Y 


( 10 ) 


where 


X-{x 1 ,X 2 >^ 3 } is the distance from the origin to the observation 
point 

Y={y^,y 2 ,y 2 ) is the distance from the origin to a source point 
and the brackets indicate that the derivative is to be evaluated at 
the retarded time t' = [t- (|X - Y|)/c]. The volume integration 
is to be performed over all space. 

It can be shown (ref. 23, e.g.) that at large distances from the 
source region (i.e., if the observation point X is many acoustic 
wavelengths from any point in the source region) and if, additionally, 
the distance between the source point and observation point is large 
compared with the dimensions of the source region (i.e., the observa- 
tion point is in the radiation field) then the spatial divergence in 
(10) becomes a time derivative and 
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V 


O 2 

, x.x. / 3 T ' . , 

p*(x,t) I — 

*«o x J 3t 




d 3 Y 


( 11 ) 


with r = |x - Y| . The argument that, at least for isothermal flows, 
T^_. is a relatively acceptable description of the acoustic sources, 
and that mean flow-acoustic interactions may be neglected as a first 
approximation was made by Lighthill and proceeds as follows. 

The term under consideration is (see fig. 2) : 


ij '0 ’ " '“i“j w ij + “ p 0 5 c 0 (p " p ofl 6 ij (12) 


t! . = (p tt + p')u.u. - e. 


If the ratio of turbulence or Reynolds stresses to viscous 
stresses is very large, that is for large Reynolds number, then 
within the region of turbulence, the term (Pq + P ' ) u j u j is very 
large compared to . Within the region of turbulence, there- 
fore, e^ is negligible. Further if the turbulent region is es- 
sentially isothermal, then the effects of heat transfer should also 
be negligible for cases where the ambient temperature is the same as 
that of the flow. Entropy spottiness within the flow should be neg- 
ligible, and the isentropic relation between pressure and density 
changes 

2 _ djp 


c = 


dp 


should hold. Consequently, within the jet (p - p Q ) % c Q (p - Pq> 


Hence the T ! . term becomes 
ij 


T Ij * Cp 0 + P ' )u i U j 

Lighthill then assumed that the density fluctuations of the flow are 


.i. . 


small compared to the mean density and hence may be neglected, so 
that finally 


T! . % p n u.u. 
ij 0 i 3 

In the region outside the flow, only acoustic disturbances need 

2 

be considered, so that p' % CqP ' , u £ u j is very small, and viscous 
effects on the acoustic waves are entirely negligible (24), except 
for propagation over very large distances. The result, therefore, 
is that Tjj is appro x imated by PqU^u. within the region of tur- 
bulence, and is zero outside the flow. On this basis, equation (9) 
can be Interpreted as a wave operator governing the propagation of 
sound whose source is given by on the right hand side. 

Equation (11) , which allows the calculation of the density 

fluctuations in the radiation field once the source term Tl. is 

ij 

known, is the usual starting point for the analysis of noise gen- 
erated by turbulent mixing, specifically jet noise. 

The present work, however, focuses its attention on noise 
generated in a turbulent combustion region. It Is certainly not 
appropriate, therefore, to neglect entropy spottiness In (12) and 
assume that d 4 - c^p’ ^ 0. Instead, in the region of combustion we 
retain the "thermal stress" term p r - c^p' as a measure of pressure 
fluctuations due to fluctuating entropy and neglect the Reynolds 
stress term. This latter assumption is made on the grounds that 
even though the turbulence intensity in a combustor may be as high 
as that in the mixing region of a jet, the flow velocities are rela- 
tively low. This would result In a relatively low level of turbulent 
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mixing noise. The quantity (p* - e^p ') , then, is retained as the 
only source term. (The viscous stress ey is still left out.) 

Still making the assumption that in the radiation field the flow 
is isentropic so that here p * ^ c^p ' , equation (11) becomes : 

p’(X,t) S'/V.t - ■f-]d 3 Y (13) 

4ttc 2 x / St 2 V V 


2 

where the source term S' = p ' - CqP 1 is a measure of the pressure 
fluctuations due to local entropy fluctuations. 

Following the procedure of Siddon for surface generated noise 
(ref. 25), equation (13) written for a time, say t*, can be multi- 
plied through by the local source pressure at Y at an earlier time 
t" - (r/cg) and time averaged: 



where the dummy variable Y* has been introduced and r’ = |x - Y f |. 
The overbars indicate a long time average. 

For flows at low Mach numbers, the retarded time differences 
can usually be neglected over the correlation volume (25) (i.e,, 
r ^ r*). Furthermore, if p’(Y) and p'(X) are individually and 
jointly stationary random variables, then their cross-correlation 
function is independent of time translations (26) : 

p^Y.t" - ^■jpCX.t') = p'(Y)p'^t + (t = t' - t") 


(15) 


T 
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Similarly, for the integrand on the right hand side of (14) 






S'(Y r ,T) (16) 


3t 


Putting Y* - Y = ri and substituting (15) and (16) into (14) 
gives : 


p’(Y,t)p'(X,t + x’) « 


P ' (Y»t) ~r S'(Y + n,t + T)d^ 


3t 


(17) 


where 


t 1 = t + ~ 


The integrand on the right-hand side of (17) is a two point, 
space-time correlation between one function and the second time 
derivative of another function, delayed in time. A theorem relat- 
ing correlations of jointly stationary random functions to their 
derivatives may be found in reference 27 , It is expressed as : 


,n .n-hn 

-H £ i (t) T5 f 2 (t + " <-« 75S £ i (t > £ 2‘ t + 

dt dt dx 


(18) 


or, in words, the cross-correlation of the time derivatives of two 
random variables may be contracted to the time derivative of the 
cross-correlation of the two variables. In equation (17) we have 
m = 0, n = 2, so (17) can be written: 


P * (Y, t)p ’ (X, t + x’) - 


3- 


p'(Y,t)S'(Y + n,t + T)dp (19) 


dx 
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If p f (Y) and S ' are uniformly coherent through the entire 
source volume, then the correlation in the integrand of the right 
hand side of (18) is independent of n» and (18) becomes: 


p'(Y,t)p'(X,t + x') !t V corr i 2P , (Y,t)S , (Y,t + T) (20) 

dx 


where V is the correlation volume: 

corr 


V = T (2i; 

corr 

For a more realistic case, where the relationship between p* 
and S’ in the source region exhibits the properties of decay in 
space and time, uniform coherence will not prevail and equation (21) 
is not an adequate representation of the correlation volume. Never- 
theless, equation (20) does suggest that as a first approximation 
we arbitrarily define the correlation volume as 


V 



p f (Y,t)S T (Y + n,t + 


x)d 3 n 


( 22 ) 


Corr p'(Y,t)S'(Y,t + x) 

With this definition, then, equation (20) becomes somewhat more 
general and will form the basis for the diagnostic phase of this in- 
vestigation. 

Subject to the assumptions made to here, (19) states that a 
cross-correlation between the pressure in a combustion noise source 
region and the pressure in the (acoustic) radiation field is propor- 
tional to the second time derivative of the cross-correlation be- 
tween that same pressure in the source region and the pressure fluc- 
tuations due to fluctuating entropy, at the same point. 
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The correlation on the right hand side of (19) may be expanded 

as: 

p'(Y,t)S'(Y,t + t) - p 1 (Y, t)jj>* (Y,t + t) - Cq P ' (Y , t + t)] 

- p’(Y,t)p'(Y,t + x) - p’(Y,t)c^p’(Y,t + t) 

(23) 

The first term on the right-hand side of (23) is simply the 
autocorrelation of the fluctuating pressure in the source region. 

The second term is a single-point cross-correlation between the 
fluctuating pressure and density in the source region. Consequently, 
if this density and pressure are in phase with each other, or if they 
are not in phase but their correlation is negligible in comparison 
with the pressure autocorrelation function, then the entire correla- 
tion on the right hand side of (19) must behave as, and therefore 
satisfy the properties of, an autocorrelation function. (If p' 
and p * are not in phase with each other, it is not likely that they 
would correlate since there is no reason to expect a systematic phase 
relationship between them.) 

For the present investigation, the two important properties are 
as follows: an autocorrelation function must be an even function 

about the t = 0 axis, i.e,, R(t) - R(~t) for all values of t; an 
autocorrelation function must have its maximum value at t = 0, i.e., 
R(0) i |R(t)| for all values of t. These two properties, then, and 
the second derivative operation on the right-hand side of (2G) , de- 
fine the shape of the correlation on the left-hand side of (20) , 
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Specifically, even order derivatives of an even function must also 
be even functions. Hence, the correlation on the left-hand side of 
(20) must be symmetric, and this symmetry must occur at a value of 
t = t 1 - r/c^, where r is the distance between the source region 
and the observation point; that is at a time delay corresponding to 
the acoustic propagation time between the source region and the ob- 
servation point. Furthermore, since the p'S* correlation should 
have a positive peak, the second derivative operation requires the 
left-hand side of (20) to have a negative peak at this axis of sym- 
metry (see fig. 3). 

Qualitatively, then, the relationship given by (20) defines the 
shape of the p'p' correlation. This can be formalised quantita- 
tively by recalling that in the frequency domain the operation cor- 
responding to time differentiation in the time domain is multipli- 
cation by frequency and a 90° phase shift. The two time derivatives 
then, produce a 180° phase shift and a multiplication by frequency 
squared. 

If the p’p* correlation in equation (23) is negligible, then 
(20) becomes: 

P* (Y, t)p'/x, t + t + — ) « p’(Y,t)p'(Y,t + t) 

\ c 0/ dr 

Taking the Fourier Transform of both sides: 

jrw/Cf, 2 

GyjjCjflOe °= -to G^-(tu) 

or 
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GyX<jw) 


-jru/c Q 

* -id e 


=3 


2 -jlr(u/c 0 )±n] 
a) e 


(24) 


where j = -/-l, and to = 2-nf (f = frequency) . 

Here, Grrr(j(u) is the cross-spec tral density between p'(Y) and 

JL A 

p' (X) , and G^p(nj) is the auto spectral density of p'(Y). 

Equation (24) is simply the definition of the transfer function 
between p'(Y), the fluctuating pressure in the source region, and 
p'(X), the fluctuating pressure at an observation point outside the 
source region: 

G-jpdw) 

h yx°“ ) 5 OjfCu) <25) 

So the transfer function between pressure in the source region 
and observation pressure has a magnitude proportional to frequency 
squared, and a phase factor given by (-rui/c^ ± it). The first part 
of the phase factor, -rw/cg, is just due to time delay and depends 
upon the distance between the two points. It produces a phase shift 
which is linearly proportional to frequency. The second part of the 
phase shift results from the second derivative operation and is in- 
dependent of frequency. These two elements of the transfer function, 
amplitude, and that part of the phase shift not due to time delay, 
uniquely determine the relationship between p'(X) and p'(Y). 

The basic approach used in the diagnostic phase of this investi- 
gation was to examine the phase and amplitude relationships between 
the fluctuating pressure in the combustion region of a turbofan en- 
gine and the fluctuating pressure at other locations within the core 







of the engine as well is in the acoustic far-field. Measured trans- 
fer functions between points within the engine which are consistent 
with the quantitative results discussed above would indicate mea- 
surable combustion noise generation (since the turbulent mixing noise 
was assumed negligible). Additionally, if the transfer function be- 
tween the pressure in the combustion region and the acoustic pressure 
in the engine far-field was also consistent with the above relation- 
ships, this would indicate that the combustion noise was, in fact, 
propagating through the engine and contributing to the overall engine 
noise. 

As will be described in the next section, the core of the en- 
gine was Instrumented with a number of pressure transducers at loca- 
tions other than in the combustion region. This was done to charac- 
terize and track the pressure signal as it propagated or convected 
through the core. To accomplish this tracking, the pair-wise time 
domain (cross correlation) and frequency domain (cross spectral 
density and coherence) characteristics were measured between points 
in the core from which noise generated by mechanisms other than com- 
bustion might be expected. 

Earlier it was mentioned that the combustion process was only 
one of a number of noise sources which could be identified with the 
core of a turbofan engine. In addition to the relatively high fre- 
quencies associated with the rotating machinery (compressor and tur- 
bine) and which can be readily removed from data by appropriate fil- 
tering procedures, there is noise generated by the flow of the gases 
through, over, and around the various internal engine surfaces and 
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passages. In general, the spectral characteristics of this noise is 
broadband and not as easily identifiable as rotating machinery noise. 
The generating mechanisms of this so called '’scrubbing” noise, or 
surface interaction noise is quite different from combustion noise, 
however. Consequently, if such noise generating mechanisms are sig- 
nificant, they should be identifiable through the appropriate trans- 
fer function as described above for combustion associated noise. 

Before beginning a description of the experimental work per- 
formed, therefore, it is appropriate to briefly describe the char- 
acteristics of these noise sources. 

Equation (11) is the solution to Lighthill's equation in the 
absence of solid surfaces. Curie (ref. 28) generalized the solution 
to Lighthill's equation to predict the acoustic radiation from a 
region of unsteady flow which contains an embedded surface. His 
solution, in the geometric and acoustic far-fields was obtained as: 


P*(X,t) 


_1_ 

4ir 




d£ + 


47TX^C r 


A [ *1 + P V. 11£ 



r 2 
3 T 

at' 


ii 


d 3 Y 


(26) 


Here, the brackets denote evaluation at the retarded time t = 
t - r/Cg as in equation (10) . The first two terms on the right hand 
side of (26) are the additional terms associated with the noise gen- 
erated by the fluid interacting with the surfaces £ , which are 


allowed to deform with surface velocity u^. The quantity is 

the local resultant stress at each point on the surface and may have 
both tangential (shear) and normal (pressure) stress components. 

The third integral is the contribution from the quadrupole noise from 
turbulence, as was shown earlier. 

Siddon (ref. 25) has shown that for cases where the volume gen- 
erated noise is weak and when the surfaces are stationary or in rigid 
steady motion (Su^/St = 0), the correlation between source pressure 
and observation point pressure is given by 

p' (X.t)p' (Y,t) « A corr ^p'(Y,t)p’(Y,t + t) (27) 

if the source pressure is uniformly coherent over the entire surface 

Here A is a correlation area analogous to the correlation 

corr 

volume in equation (21) . If the source pressure is not uniformly 

coherent, then A may be generalized analogously to V in 

corr corr 

equation (22), with similar restrictions. 

It is clear from equation (27) that the transfer function be- 
tween source pressure p’(Y) and observation point pressure p'(X) 
can be obtained similarly to equation (24) as 

S^rCjw) -jt(r<u/c n ) + 0f/2)] 

"a 010 5 

For surface generated noise, then, the magnitude of the trans- 
fer function between source pressure and observation point pressure 
is linearly proportional to frequency and the part of the phase shift 
not due to time delay is ±90°, independent of frequency. In the time 


domain, the cross-correlation is the first derivative of an even 
function, and so it will be an odd function, with the axis of anti- 
symmetry being at a value of t = r/cg (see fig. 3). 

We see, then, that the phase and amplitude relationships which 
prevail between the pressure in a source region and in a non-source 
region are uniquely determined by the acoustic generating mechanism. 
The transfer functions, which can be written in terms of the phase 
and amplitude between the various pressures can then be used to 
identify this mechanism. 

Finally there exists the possibility that neither pressure 
transducer of a pair within the engine which is being correlated is 
in a source region. In this case the transfer function between the 
two should have a constant magnitude (i.e., flat and independent of 
frequency) and a phase factor due to time delay only. The relative 
phase between the two which is not due to time delay should be zero. 

These relatively straightforward notions regarding transfer 
functions between pressures within the engine core and/or outside 
the engine (i.e., the acoustic pressure) form the basis for the 
diagnostics to be described. 


3 . DESCRIPTION OF EXPERIMENT 


3 . 1 Engine 

The engine which was the object of the current investigation 
was an AVCO-Lycoming YF-102 turbofan engine. The core, or gas 
generating section of the engine was originally designed for and 
used to power gas turbine driven helicopters. The fan stage was 
added by Lycoming to convert the engine into a turbofan. This en- 
gine will eventually be the propulsion system for the NASA Quiet 
Short-Haul Research Aircraft (QSRA) , with the airframe designed and 
built by The Boeing Company. 

A cutaway illustration of this engine in the unmounted con- 
figuration is shown in figure 4. The engine has a relatively high 
by-pass ratio of six to one and a rated thrust of 33 kN. The gas 
generating section consists of an eight stage compressor, seven 
axial and one centrifugal, driven by the first two (high pressure) 
stages of the turbine with a reverse-flow annular combustor. The 
fan and supercharger are driven through a 2.3 to 1 speed reduction 
gear by the second two stages (low pressure) of the turbine. The 
fan is 1 meter in diameter, the primary nozzle exit area is 0.136 
square meters (42 cm diam) and the fan nozzle exit area is 0.361 
square meter. The engine was tested with a bellmouth inlet to pro- 
vide for smooth inflow conditions. A photograph of the engine 
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mounted on the test stand is shown in figure 5. 

It should be noted that the reverse flow-combustor in this en- 
gine is an unusual arrangement for most aircraft power plants , es- 
pecially those in commercial service. Essentially the design repre- 
sents a compromise made to minimize engine length. As such, it is 
more appropriate for helicopters and small business jet aircraft. 

It is possible that the characteristics of combustion noise gen- 
erated by the annular or can-annular combustors on commercial turbo- 
fan engines without reverse flow are different than those of the 
reverse-flow combustor. This may be particularly true of the spec- 
tral content, which is likely to be sensitive to combustor geometry. 

The final point to be made about the YF-102 engine concerns its 
moderately high by-pass ratio of six to one, which permits a rela- 
tively low primary jet exhaust velocity of approximately 315 
meters/sec at takeoff thrust. This compares with approximately 
550 meters/sec primary jet velocity for the Pratt & Whitney JT8D 
turbofan which has a bypass ratio of 1.05 to one. Because of the 
much lower jet velocities for the YF-102 engine, the relative combus- 
tion noise contribution to overall engine noise may be greater than 
for an engine like the JT8D, which powers approximately 65 percent 
of commercial jet aircraft in this country. Consequently, for these 
reasons, caution should be exercised in generalizing the results to 
be presented later to the lower bypass ratio turbofan engines with 
more conventional combustor configurations. 
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3.2 Internal Pressure Instrumentation 

The fluctuating pressure measurements within the engine core 
were made at seven different locations* as shown in figure 6. Their 
number and locations were: two just downstream of the compressor 

exit about 2 cm apart; one at the combustor entrance; two within 
the annular combustor itself, both at the same axial location but 
separated 90° circumferentially; and two within the primary nozzle, 
one just downstream of the turbine at the primary nozzle entrance, 
and one close to the nozzle exit plane* 

This thesis, however, will report the results obtained from 
the two combustor probes and the two tailpipe probes only. The 
three other pressure probes were used to obtain data in a parallel 
program at the Lewis Research Center and will not be discussed here. 

The actual transducers used were conventional 0,635 cm condenser 
microphones with pressure response cartridges. To avoid direct ex- 
posure of the microphones to the high temperature combustion gases 
within the core, they were mounted outside the engine and the fluc- 
tuating static pressure in the engine core was communicated to the 
transducers by means of a "semi-infinite" acoustic waveguide tube. 

As described below, this waveguide tube consists, essentially, of a 
very long extension of narrow tubing. The long length is required 
to make the tube essentially anechoic to eliminate longitudinal 
probe resonances (i.e., standing waves). 

A drawing of a typical probe is shown in figure 7. The micro- 
phone was flush mounted in the acoustic waveguide through a support- 
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ing block and housed in a pressure chamber. Attached to the block 
were a 5/8-cm diameter sensing tube on one end and a coil of tubing 
of the same diameter, 30 meters long on the other. The diameter of 
the coil was approximately 30 cm, so the ratio of the tube radius to 
coil bend radius was approximately 50 to 1. This was considered suf- 
ficient to minimize the possibility of reflections from the bends in 
the tube. 

The lowest frequency at which transverse resonances (i.e., 
radial or circumferential) could be expected was estimated from 
(ref. 29) 

„ 1.84 ire 

f = 5 

Where c is the speed of sound in the tube, and D is the tube 
diameter. For c = 340 m/sec, and D = 5/8 cm, then f = 31.8 kHz. 
This is well beyond the frequency range of interest, as will be 
seen. 

The sensing tube of each probe was flush mounted as a static 
pressure tap at each of the various measuring locations within the 
engine core. A regulated nitrogen purge flow was maintained in each 
sensing line to protect the microphone from hot core gases. Since 
the signals from various pairs of probes were to be cross- 
correlated, the regulated purge flow system was separate for each 
probe* This eliminated the possibility of any common valve noise 
from the regulator contributing to the correlations* Static pres- 
sure was balanced across the microphone by means of a small vent hole 
connecting the pressure chamber and sensing line. Figure 8 is a 
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schematic of a typical core pressure probe installation. 

An inescapable characteristic of such probes is the frequency 
response roll-off with increasing frequency due to viscous attenua- 
tion in the tube between the sensing opening and the microphone. 
Prior to the tests, the frequency response of the probes were deter- 
mined. The amplitude and phase response were determined by compar- 
ison with a microphone identical to the one used in the probe using 
a symmetric placement with respect to the axis of a loudspeaker, in 
an anechoic room. The input to the loudspeaker was a signal from a 
white noise generator low pass filtered at 10 kHz. 

The frequency response for a typical probe is shown in fig- 
ure 9. The results indicated that the amplitude response of the 
probes was flat within ±2 dB from 50 to 1500 Hz. Similarly, the 
phase response of the probes was flat within about 5° up to 1500 Hz 
after accounting for the phase lag associated with the length of 
the sensing tube of each particular probe. Between 1500 Hz and 
about 3500 Hz the response was generally flat within ±4 dB and 
±10°. Beyond 3500 Hz the response curves showed large distortions, 
and it was not clear whether these distortions were due to limita- 
tions in the calibration facility or were intrinsically related to 
probe construction and geometry. However, even at 1500 Hz, with the 
exception of the high frequency tones generated by the rotating 
machinery, the core noise associated with this engine was confined 
to a frequency range well within the acceptable response region of 
the probes. 
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A facility for probe calibration at the elevated temperatures 
prevailing within the core of the engine (between 810 K and 1375 K, 
depending upon operating conditions) was not available. Consequently, 
the actual response of the probes during testing cannot be docu- 
mented. However, it is felt that the presence of the nitrogen 
purge, the primary purpose of which was to prevent a blow-by of hot 
gases across the microphone during engine start-up and operating 
transients, served to significantly cool the sensing tube. The 
temperature within the tube, just upstream of the microphone (see 
figs. 7 and 8), was monitored with a thermocouple. During engine 
operation the temperature did not exceed ambient by more than about 
5 K, 

The purging nitrogen flowing through the transmission tube can 
generate pseudosound as it passes across the face of the microphone. 
This pseudosound pressure level was measured directly with the en- 
gine off and found to be a minimum of 20 dB below the fluctuating 
pressure levels in the core at the highest flow rates required dur- 
ing testing. The nitrogen purge flow rate was monitored (see fig. 8) 
during engine operation for the pressure level comparison between en- 
gine on and engine off. 

The probe sensing tubes and anechoic extensions were required 
to cross the fan flow stream when mounted on the engine. Impinge- 
ment of the fan flow on the tubes, then, would induce noise within 
the tubes themselves, as well as radiate noise to the far-field. To 
minimize this flow Induced noise, airfoil shaped struts attached to 
the cowl of the core engine and surrounding the tubes were used. A 
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photograph of the engine with the probes and struts in place is 
shown in figure 10. 

3.3 Test Arena and Far-Field Microphone 

Array; Test Conditions 

All tests were conducted at an outdoor test facility with a 
hard-surface (concrete) ground plane, at the NASA Lewis Research 
Center. The engine was suspended from the test stand with its 
centerline 2.9 meters above the ground plane. 

The far-field microphone array consisted of sixteen 1.27 cm 
diameter condenser microphones placed on circular arcs centered on 
the pivot plane of the test stand (i.e., the plane at which the 
engine is suspended from the test stand) . The microphones were 
spaced at intervals of 10° from 10° to 160°, relative to the en- 
gine inlet axis. To minimise problems associated with ground plane 
reflections, the microphones were mounted at ground level. Addi- 
tionally, all microphones were fitted with wind screens. Because 
of the presence of small utility junction boxes on the ground, the 
radius of the rearward arc (29.9 m) was 0.6 meter less than the 
radius of the forward arc (30 ,5m). Figure 11 is a schematic il- 
lustration of the test arena. 

The simultaneous internal (i.e., core) fluctuating pressure and 
far-field acoustic measurements were made at eight different fan 
speeds at approximately equal intervals between 30 and 95 percent 
of maximum speed (7600 rpm) . The corresponding range of combustor 
temperatures and primary jet exhaust velocities were from 810 K and 
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98 m/sec to 1375 K, and 314 in/ sec. The ambient temperature during 
testing ranged between 0° and 13° C, and the relative humidity ranged 
between 50 and 70 percent. Data were not taken during periods of 
precipitation or when the wind velocity exceeded about 10 km/hr 
(5 knots) . 

3.4 Data Acquisition 

Prior to operation, all microphones were calibrated with a 
pistonphone which provided a standard sound pressure level of 
124 dB (referenced to 2x10 ^ Pa) at 250 Hz. During testing, sig- 
nals from all transducers were transmitted to a remote control room 
via shielded cable for amplification, signal conditioning, and re- 
cording . 

The signals from the far-field microphones and the core fluc- 
tuating pressure probes were FM recorded on magnetic tape for later 
processing. The tape recorders were equipped with IR.IG standard, 
intermediate band, group II electronics, which provide a flat fre- 
quency response to 20 kHz on both record and reproduce at a tape 
speed of 60 in. /sec. The dynamic range available was approximately 
48 dB. Two fourteen channel tape recorders were used with all seven 
internal probe signals recorded on each machine, with the 40° 
through 100“ far-field signals recorded on one, and the 110° 
through 160° far-field signals on the other. Additionally, one 
channel on each machine was used for annotation. The 10°, 20°, and 
30° microphone signals wi re not recorded simultaneously with the 
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remainder of the data and so could not be used for correlation pur- j 

poses , 

All internal fluctuating pressure signals and far-field acoustic j 

signals were recorded in 120 second record lengths. This was init- ! 

ially considered sufficient to obtain high statistical confidence j 

; i 

in the data. However, in some cases it provided only marginal ■■ j 

j 

confidence levels, and these will be noted later, u 

Engine operating conditions as well as probe temperature, pres- 

l i 

sure, and nitrogen purge flow rate were monitored and recorded on j 

the NASA-Lewis CADDE system (Central Automatic Digital Data En- 
coder) . On-line data processing and analysis was limited to 1/3 
octave-band spectral analysis of the various pressure signals for 
"quick-look” purposes. Obviously anomalous data were discarded. 

The Lewis Research Center is immediately adjacent to Cleveland 
Hopkins International Airport, and occasional aircraft flyovers 

j 

would require re-running a data point. jv 

3.5 Data Processing 

Virtually all of the results presented herein were obtained by 
off-line data processing of the taped signals on a two-channel, FFT 
(Fast Fourier Transform) digital signal processor, with built-in 
analog to digital converters and 120 dB/octave anti-aliasing fil- 
ters. The processor was capable of direct computation of up to a 
4096 sample ensemble average of a 1024 point forward or inverse 

Fourier transform to yield either single or double channel time _ 
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domain (auto or cross-correlation) or frequency domain (auto or 
cross-spectra, phase, transfer function, and ordinary coherence) 
information. 

The processor also had a unique feature which permitted manual 
editing of stored input data prior to processing. For example, the 
input memory could have its contents rotated relative to a fixed 
memory address register. If the stored data represents a time 
series, then rotation represents a time translation. A Fourier 
transform could then be performed on the rotated data. This very 
useful feature enabled a previously computed cross-correlation func- 
tion to be rotated in memory for time delay removal. A Fourier 
transform of the translated cross-correlation then produces the 
relative phase shift between the two correlated signals unobscured 
by that part of the phase shift due to time delay. 

There was a slight limitation in the selection of anti- 
aliasing frequencies imposed by the design of the processor. To 
understand this limitation, a brief description of the frequency 
selector settings of the processor is necessary. The analyzer had 
available 57 front panel selectable frequency ranges, from 10 Hz 
full scale to 150 kHz full scale in integer steps of 1 through 15, 
with multipliers of 10, 100, 1000, and 10,000. These 57 frequency 
settings represented the 0 Hz (nominal) to full scale frequency 
analysis ranges available to the user. However, the anti-aliasing 
filters could only be activated at cutoff frequencies equal to 
exactly 80 percent of any one of the full scale analysis ranges (not 
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necessarily the range being used) . 

For example, if it was desired to analyze data in a frequency 
range of 0 to 400 Hz, then the anti-aliasing frequency could be 
selected at precisely 400 Hz since 400 Hz is 80 percent of 500 Hz, 
and 500 Hz was a selectable frequency range. However, if it was de- 
sired to analyze data in a frequency range of, say, 0 to 1000 Hz, 
it was not possible to select an anti-aliasing frequency of ex- 
actly 1000 Hz since this would have required a frequency range 
setting of 1250 Hz to be available, which it wasn't. Hence, in 
this case, an anti-aliasing frequency of 1040 Hz would be required 
since 1300 Hz as a frequency selection was available. Consequently 
some of the data contains small aliasing errors due to this char- 
acteristic of the analyzer. This is considered to be a relatively 
small compromise in return for the avoidance of using external fil- 
tering with lower roll-off rates, and the possible introduction of 
small phase shifts due to mismatched filters. At any rate, in all 
the data to be presented, both the frequency range of analysis as 
well as the anti-aliasing frequency selected will be specified. 

Another limitation in the processor, this one somewhat more 
serious than above, concerns the computation of coherence functions 
between signal pairs. In the usual analysis of signal pairs, the 
similarities or mutual properties between the two signals do not 
coincide with one another until some time delay has elapsed. If 
this time delay is on the order of the sampling period which is 
being used for analysis, this will lead to a low estimate of the 
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true coherence function. In a fixed point analysis, however, the 
sampling period is related directly to the frequency analysis range 
and is not independently controllable. 

The following illustrative example will serve to clarify this 
point. The signal processor used to analyze the data in this in- 
vestigation digitizes each time record into 1024 words at a sampling 
rate of 2.048 times the highest frequency selected (i.e., the max- 
imum frequency in the range selected. Note that this is slightly 
in excess of the Nyquist sampling criterion) . Hence 

f , = 2.048 f 

sample m 

where f is the maximum frequency to be resolved, 
m 

The memory period, or sample record length, therefore is: 

T _ Number of words/record 
Number of words/ sec 


= 1024/(2.048 f ) = 500/ f sec/record 
m m 

If, for example, frequency information to 10,000 Hz is re- 
quired, then f — 10,000, and T = 500/10,000 = 0.05 seconds. If 
m 

the natural time delay between the two signals being analyzed is on 
the order of 0.05 second or greater, than a significant bias error 
is introduced into the measured, or estimated coherence functions, 
with the estimate being low (30) . 

The required procedure to eliminate this bias error is to in- 
troduce a pre-computation delay (i.e., sample and hold) In the 
signal which occurs earlier in time, thus waiting for the second 
signal to "catch up" before Implementing any computations. The 


processor had such a precomputation delay feature when computing 
cross-correlations and cross spectra, but not when computing co- 
herence functions. This did not present any problems when com- 
puting coherence functions between pressures within the engine, 
since the natural time delays were relatively short in comparison 
with the processor memory periods. In any event, the in engine 
results were used for diagnostic purposes only. 

Between the engine and the far-field, where separations of 
approximately 30 meters and 90 msec prevails, however, the time 
delays were a significant fraction of the processor memory period 
for frequency analysis ranges above about 500 Hz. Furthermore, 
it was between the engine and the far-field that the numerical 
value of the coherence was required for computing far-field com- 
bustor spectra, directivity and power levels. For all cases where 
the numerical value of the coherence was used for computations the 
bias correction due to time delay was calculated, and its effect 
on the resulting computation noted. 

Finally, it should be pointed out that the somewhat less than 
flat response of the waveguide probes do not affect the quantitative 
results to be shown in Chapter 5. The normalized nature of the co- 
herence function is such as to cancel the effect of probe response 
(as long as the response is linear) , All measurement errors , there- 
fore, are determined by the microphone response themselves, which 
generally are accurate to better than +1 'IB. 


4. ^IDENTIFICATION OF COMBUSTION NOISE 



4.1 Correlations Between Internal Pressure 

This chapter discusses the results of the diagnostics con- 
ducted in the present investigation to identify combustion noise in 
the YF-102 turbofan engine. We begin by noting that when conduct- 
ing fluctuating measurements in a turbulent fluid for purposes of 
obtaining direct acoustic information, the investigator is nor- 
mally confronted with a significant problem in interpreting his 
data. There is no straightforward way of distinguishing between 
that part of the pressure signal which is acoustic information, and 
that part which is dt^ to hydrodynamic pressure fluctuations arising 
from local or convecting turbulence. 

In general, however, the convection speeds of turbulence 
within an engine are much less than the local speed of sound. This 
leads to the idea of cross-correlating the pressures at two sepa- 
rate points within the engine and examining the time delay. If 
the pressure signal is entirely due to convecting turbulence, then 
the cross-correlation will be a positively peaked, even function 
with respect to a time delay corresponding to the convection speed. 
If the pressure signal is due entirely to a propagating acoustic- 
wave, then the cross-correlation will have similar properties with 
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respect to a time delay corresponding to the acoustic propagation 
speed (augmented or diminished by the speed of the flow, depending 
on direction) . 

If the pressure signal is due to both acoustic waves and con- 
vecting turbulence then the cross— correlation should show two posi- 
tive peaks, with each peak occurring at the appropriate time delay, 
as described above. For this case, digital computational tech- 
niques have been developed which utilize the symmetry properties 
of correlation functions to separate the acoustic contributions to 
the correlation from the contribution due to convecting turbulence 
(ref. 31). This may then be Fourier transformed to obtain the 
acoustic spectrum. It is not necessary to remove the time delay, 
since time delay produces phase shifts only, and does not alter 
the spectral shape. (It should be noted that this technique re- 
quires the acoustic information to be uncoupled from the turbu- 
lence. If the two are not independent, then their individual cross- 
correlations are not separable in this manner. In aerodynamic noise 
problems, of course, the acoustics are very often coupled to the 
turbulence . ) 

With the objective of characterizing the internal fluctuating 
pressures in the tailpipe of the YF-102, these two pressure signals 
(see fig. 6) were cross— correlated. The normalized result is shown 
in figure 12^. The correlation was obtained at a fan speed of 

^All correlation data to be shown in this chapter will have the 
following notation shown on the figures: The anti-aliasing filter 

cutoff frequency, f c ; the number of disjoint (i.e., independent) 
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43 percent of maximum, which is near flight idle conditions , but 
below approach conditions. Both signals were low-pass filtered at 
1600 Hz with the anti-aliasing filters of the digital signal proc- 
essor (DSP) . This cutoff frequency was chosen as the approximate 
upper limit of the acceptable response range of the probes. As will 
be shown shortly, this is not a significant limitation since the 
combustion noise for this engine is confined to low frequencies. 

The correlation was normalized with respect to the zero time 
delay values of the individual auto-correlations (i.e., the in- 
dividual rms values within the 1600 Hz passband) . The signals 
correlate rather strongly at a positive delay time of 0.54 msec. 

The separation between the two tailpipe probes was 0,34 m. Conse- 
quently, considering the gas temperature within the nozzle (^670 K) 
and the flow velocity (^90 m/sec) for this case, the calculated time 
for an acoustic wave to travel between the two probes is about 
0.56 msec. At the speed of the flow, the corresponding convection 
time would be about 3.8 msec. 

This close agreement between the measured delay time and the 
calculated acoustic delay time must be considered fortuitous. Al- 
though, as indicated earlier in Chapter 3, the temperatures in the 
probe sensing tubes were monitored and were within about 10 K of 


samples used for ensemble averaging, N, and the digitizing time in- 
crement, At. Additionally, frequency domain data will have the fol- 
lowing notation on the data-f igures s Upper frequency of analysis 
range, f m ; anti-aliasing frequency, f c ; number of disjoint samples 
for the ensemble average, N; resolution bandwidth, Af. Correlation 
data contain 1024 points. Frequency domain data contain 512 ampli- 
tude or 512 phase points. 


ambient, some temperature gradients near the open end of the tubes 
must exist. These temperature gradients could introduce some small 
relative phase shifts between probes. Any such phase shifts, in 
turn, would show up as time delay between the signals , For the 
small measured time delays between signals within the engine, then, 
conclusions regarding the actual speed of the signals are some- 
what speculative. Nevertheless, this correlation appears to be 
evenly symmetric with a positive peak about the acoustic delay 
time, and the preliminary evidence suggests, that the two pressure 
probes are detecting an acoustic signal. 

Before proceeding with additional results, two points should 
be made regarding the above conclusion. First, under certain con- 
ditions a cross-correlation which has a single peak at a delay 
time which corresponds to the speed of some identifiable phenom- 
enon, and which meets the symmetry conditions discussed above, 
does not ensure that that phenomenon alone is, in fact, being 
measured. The width of the main lobe in the correlation of a 
broadband random signal is approximately equal to 1/B W , where B w 
is the bandwidth of the signal (26) . Clearly then, unless the dif- 
ference in time delays being investigated is substantially greater 
than 1/B w , the two separate peaks may not be distinguishable, due 
to finite signal bandwidth. In the present case the two time de- 
lays being investigated differ by about 3.2 msec, which is approx- 
imately five times the reciprocal of the signal bandwidth, A 
second peak should be clearly distinguishable If It were present. 
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More important is the fact that if the upstream acoustic signal 
was contaminated by turbulence which decayed before reaching the 
downstream probe (the signal at which may or may not be contami- 
nated) , the resulting correlation would still show a single peak. 

The conclusion, then, that the pressure signal within the nozzle is 
a purely acoustic signal can only be tentative. The issue of tur- 
bulence contamination on the various pressure signals will be dis- 
cussed later in this chapter. 

Returning to the results, if it is assumed that the nozzle 
pressure is related to an acoustic wave traveling in the downstream 
direction, then it is natural to examine the pressures upstream of 
this region to determine the origin of this acoustic information. 

The next upstream probe location is the combustor. The two com- 
bustor probes, 90° apart circumferentially but at the same axial 
location, were flush mounted in the combustor liner downstream of 
the igniter in the region near where the combustion process is 
completed (see fig. 6). 

The normalized cross correlations between the pressure at the 
in-line combustor probe location and the pressure at each of the 
two tailpipe positions for the same engine conditions as fig- 
ure 12, are shown in figures 13(a) and (b). These correlations 
are seen to be significantly different from that of figure 11. 

They contain both positive and negative multiple peaks and are more 
complexly shaped than the cross-correlation shown in figure 12* 

They clearly do not have, for example, the symmetry and maxima 
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properties associated with pure propagation or convection (i.e., 
just time delay) • 

This result can he clarified considerably by examining the 
frequency domain information between these pairs of pressure sig- 
nals. 


4 . 2 Pressure Coherence 


The tool which is used here to aid in interpreting these corre- 
lations is the ordinary coherence function. The ordinary coherence 
function is completely analogous to the cross-correlation function, 
with the information being presented in the frequency domain. It 
is essentially a normalized cross-spectrum and is defined for random 
signals as (26) : 
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where | G^ b (ju) | is the square of the ensemble averaged cross- 
spectral density between a and b; and G (to) and G bb (iu) are 
the averaged autospectral densities at a and b, respectively. 

The coherence function must have a value between zero and one, 
with high (low) coherence at a particular frequency, f, meaning 
high (low) correlation at that frequency. Of course, here time 
domain information has been traded for frequency domain information. 
Subject to certain limitations, the coherence function may also be 
used to compute quantitative contributions between a source region 
and a field point. These limitations, and the results of such com- 
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putations will be presented in the next chapter. 

The measured coherence function between the pressure at the in- 
line combustor station and the pressure at the upstream nozzle sta- 
tion is shown in figure 14. There are three very distinct and 
separata regions of coherence: one between zero and 250 Hz; another 

between about 350 Hz and 650 Hz; and a third between about 750 Hz 
and 950 Hz. The issue of the numerical value of the coherence will 
be dealt with in the next chapter. Here we are concerned only with 
the fact that there are three separate spectral regions where the 
fluctuating pressure in the combustor is coherent with the pressure 
in the upstream region of the tailpipe. 

Particular attention should be paid to the second of these 
three regions, between about 350 Hz and 650 Hz. This is precisely 
the region where numerous investigators (refs. 13 to 19, for 
example) have found the spectral range of core noise to lie. This 
frequency range has been the demonstrated concern among most inves- 
tigators of core noise. In fact, at least one of the more widely 
used prediction schemes considers the spectral peak of core noise 
to occur at a fixed frequency of 400 Hz, independent of engine 
geometry or operating conditions (32) . 

In figure 15, however, which shows the measured coherence be- 
tween the fluctuating pressure in the combustor and the pressure at 
the downstream tailpipe station, this central region of coherence 
has diminished significantly, while the low frequency region of co- 
herence is still clearly present. The highest frequency region of 




coherence has vanished completely. 

Similarly, figure 16 shows the measured coherence function be- 
tween the fluctuating pressure in the combustor and the far-field 
acoustic pressure, at the 120° microphone. This microphone loca- 
tion was selected since data taken by other investigators using 
more conventional methods appear to indicate that the angle of peak 
core noise occurs near 120° (see ref. 32, for example) . This would 
improve the likelihood of obtaining a relatively strong coherent 
relationship between the two measuring stations if the core noise 
had a significant combustor related component. As can be seen from 
the figure, the coherence between the combustor pressure and the 
far-field acoustic pressure is restricted entirely to the same low 
frequency regions of coherence (0-250 Hz) which prevailed between 
pressures within the engine. 

This result indicates, then, that insofar as combustor asso- 
ciated noise is of interest, the regions of fluctuating pressure 
coherence above about 350 Hz which are seen in figures 14 and 15 
are not associated with sound, or with a linear sound producing or 
propagating mechanism. Alternatively, as shown in figure 16, if 
these pressure coherence regions are associated with sound, the en- 
ergy contained in those spectral regions is not radiating outside 
the engine. Consequently, for purposes of examining combustor re- 
lated far-field noise in this engine, attention can be limited en- 
tirely to the frequency range below about 250 Hz, which clearly is 
radiating to the far-field. 
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4.3 Phase and Amplitude Relations 

It is now of interest to reexamine the cross-correlation shown 
in figure 13 (a) , allowing the coherence result shown in figure 16 
to serve as a guide for filtering the data. The unnormalized fil- 
tered cross-correlations between the combustor fluctuating pressure 
and the nozzle fluctuating pressures (i.e., the filtered versions 
of figs. 13(a) and (b) , respectively) are shown in figure 17. All 
signals have been low-pass filtered at 240 Hz, the frequency at 
which figure 16 indicates there is no combustor related far-field 
sound. The correlations were not normalized since their magnitudes 
are not important to the present discussion. 

As can be seen, these cross-correlations are much better char- 
acterized than the broadband versions shown in figure 13. They 
have clearly defined negative peaks evenly symmetric about delay 

times of about 5.1 msec in figure 17(a) and 5.7 msec in figure 

2 

17(b) . This, of course, is one of the characteristics suggested 
in the discussion In Chapter 2 which is expected to prevail be- 
tween the pressure in a combustion noise source region and the 
acoustic pressure outside such a source region. The difference in 
these two delay times of 0.6 msec corresponds very closely to the 
measured acoustic delay time of figure 12, as would also be ex- 
pected. 

2 

This kind of precision is obviously not obtainable from the 
figures, which were drawn on an x-y plotter driven by the DSP. 
The result was obtained by direct digital readout of the output 
memory of the DSP. 
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No conclusions regarding the absolute delay times, however, 
can be made. The flow path between the combustor probe and the up- 
stream tailpipe probe is exceedingly complex as can be seen from 
figures 4 and 6. Furthermore, neither the temperature distribution 
nor flow velocities throughout this region were known. Conse- 
quently, whether the observed 5 msec delay time corresponds to an 
acoustic propagation time or a turbulent convection time is not at 
this point, obvious. In either case, however, because of the nega- 
tive peaks, these correlation functions cannot be associated with 
pure time delay. This reasoning is based on the notion that if only 

pure time delay prevails between two signals then the cross- 

correlation should be shaped like the autocorrelation, which must 
have a positive peak, translated in time. 

The nature of the relationship between the combustor and tail- 
pipe pressures in this low frequency range of combustor associated 
far-field noise can be further explored by examining the relative 
phase difference between the two signals. The phase shift between 
combustor and upstream tailpipe pressures, as measured directly 
with the DSP, is shown in figure 18, for the signals low pass fil- 
tered at 240 Hz. The plot shows a phase shift between the two sig- 

nals which is linear with frequency, and which Is characteristic of 
the phase relationship between two signals with time delay between 
them. 

The shape of a cross-correlation function is entirely deter- 
mined by the amplitude relationship between the two signals and by 
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that part of the phase relationship which is not due to time delay. 
Time delay causes a phase shift which is linear with frequency and 
which serves only to translate the correlation function. What is 
desired, in this case, is the phase difference unobscured by time 
delay . 

This zero time delay phase shift corresponding to figure 18 
was obtained by removing the time delay associated with the fil- 
tered correlation function between combustor and upstream tailpipe 
pressure (fig. 17(a)) via the editing feature of the DSP described 
in Chapter 3 and Fourier transforming the translated correlation. 

The real and imaginary parts were then combined in polar form to 
produce the desired phase information (as well as the amplitude of 
the cross spectrum). However, there are some subtle, but impor- 
tant, features associated with such an operation. 

The correlations which have been presented to this point have 
been computed by low pass filtering the signals at some appropriate 
frequency but sampling at a much higher frequency. For example, 
figure 18 was obtained by filtering the data at 240 Hz but selecting 
an analysis range of f m = 7000 Hz. Since the memory period for 
this DSP is 500/f m the equivalent time scale for the resulting 
cross-correlation is 500/7000 Hz = 71.4 msec. Hence each of the 
1024 points in the correlation represents a time interval of 0.0698 
msec. Sampling the data at a frequency much higher than the filter 
frequency, as in this case, has the effect of expanding the time 
scale and presenting the correlation function in great detail, so 
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that any symmetry properties can be easily observed and the time 
delay precisely located. 

However, since the DSP uses a fixed number of points (1024) in 
performing an FFT on the data, if this expanded correlation is 
Fourier transformed, the 240 Hz information would be compressed in 
the frequency domain. This compression factor would be, in this 
case, 240 Hz/7000 Hz = 0.0342, so only about 3.4 percent of the re- 
sulting information would be useable, or 17 phase points (and 
17 amplitude points). This makes the data difficult to interpret 
and also increases the effective resolution bandwidth with the pos- 
sible introduction of bias errors. Consequently, the correlation 
function must be recomputed at a lower sampling rate, and then 
edited (i.e., rotated in memory) to remove the time delay. 

The result of this computation at a lower sampling rate (but 
still low-pass filtered at 240 Hz) for the combustor to upstream 
correlation function is shown in figure 19. Comparing this with 
figure 17 (a) , it is clear that now the time domain information has 
been compressed with some loss of detail. However, the symmetry 
properties which suggested that it is the time delay to the nega- 
tive peak which is relevant, as well as the value of the time delay 
was obtained from figure 17 (a) . This compressed correlation was 
computed using an analysis range of 1000 Hz, so now the compression 
factor is 240/1000 = 0.24. Now if the time delay to the negative 
peak is removed and the result Fourier transformed, some 24 percent 
of the resulting frequency domain information will be useable. 


The phase shift obtained by translating the negative peak of 
the correlation in figure 19 to zero time delay, and Fourier trans- 
forming is shown in figure 20. It is seen here, that with time de- 
lay removed, the phase shift between combustor pressure and up- 
stream tailpipe pressure is essentially 180°, independent of 
frequency, for the 240 Hz low pass band of Interest. 

The relative amplitude relationship between the combustor 
pressure and the upstream tailpipe pressure can be obtained via the 
transfer function. The relative amplitude of the measured transfer 
function (not the cross spectrum) between these two pressure, com- 
puted directly by the DSP is shown in figure 21. In the present 
context, this transfer function is defined as the ratio of the 
cros 5-spectrum between the two pressures, normalized with respect 
to the auto-spectrum of the upstream pressure. That is, it is 
precisely the transfer function of equation (25) . Comparing the 
result with the dashed line sketched on the figure, it is seen that 
the transfer function between combustor pressure and upstream tail- 
pipe pressure varies approximately as the square of the frequency 
(6 dB/octave), over most of the range of interest. 

A least squares linear curve fit of the transfer function 
shown in figure 21 between 10 and 150 Hz yielded an actual fre- 
quency exponent of 1.7 (^5.2 dB/octave). There are several pos- 
sible reasons why an exponent of precisely 2.0 was not obtained, 
but the most likely is the presence of some "scrubbing" noise 
within this frequency range. As discussed earlier, scrubbing, or 
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surface interaction, noise would be characterized by a transfer 
function proportional to frequency to the first power. A small 
amount of scrubbing noise, then, would serve to reduce the combus- 
tion noise frequency exponent somewhat. In fact, all of the 
transfer functions to be shown between the combustor pressure and 
downstream pressures, including far-field, have actual frequency 
exponents somewhat less than 2.0. However, all of the correspond- 
ing cross-correlations and time delay removed phase shifts to be 
shown between combustor and tailpipe pressures and between com- 
bustor and far-field pressures exhibit the symmetry properties and 
180° phase shift, respectively, as discussed earlier for combus- 
tion noise. 

Beyond about 150 Hz, figure 21 shows a sharp decrease in the 
transfer function. Similar behavior in all the transfer functions 
measured with respect to the combustor pressure was observed, as 
will be seen shortly. This would indicate an attenuation mech- 
anism downstream of the combustor strongly dependent on frequency, 
and beginning near 150 Hz. Such a result was reported in refer- 
ence 12 at frequencies only slightly higher than the 150 Hz shown 
here, and was attributed to turbine attenuation. It is concluded, 
therefore, that although an exact frequency squared relationship 
is not observed in figure 21, there is adequate evidence to support 
the contention that the dominant behavior is consistent with equa- 
tion (24) for the combustor as an acoustic source region. 

Similar measurements were made to obtain the phase and ampli- 
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tude relationships between the combustor pressure and the down- 
stream tailpipe pressure. The broadband cross-correlation was 
shown in figure 13(b) . The filtered cross-correlation over the 
240 Hz low pass band of interest was shown in figure 17(b). The 
phase, before time delay removal, as obtained by direct measure- 
ment on the DSP is shown in figure 22(a). Again we see a phase 
shift linearly proportional to frequency. The phase shift after 
time delay removal, obtained by Fourier transforming the com- 
pressed correlation function (not shown) is shown in figure 22(b), 
and here we also see a phase shift of 180°, essentially indepen- 
dent of frequency, over most of the frequency range of interest. 
Finally, the magnitude of transfer function between the combustor 
pressure and the downstream tailpipe pressure was also measured 
with the DSP and this is shown in figure 23. The result, once 
more, is a relative amplitude difference approximately proportional 
to frequency squared. 

At this point, it is worth summarizing the results obtained 
so far. It has been shown by coherence measurements that the 
fluctuating pressure in the combustor correlated with the pressures 
in the tailpipe over several distinct spectral ranges. However, 
the coherence between the fluctuating combustor pressure and the 
far-field acoustic pressure is restricted entirely to a single re- 
gion below about 250 Hz, and peaking near 125 Hz. That is, only 
the low frequency range of the fluctuating combustor pressure is 
related to far field noise, the subject of interest. 
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Using this result as a guide, filtered cross-correlation func- 
tions between the combustor pressure and the tailpipe pressures have 
been measured. These filtered correlation functions exhibit the 
property of being evenly symmetric about negative peaks and hence 
cannot be associated with pure time delay (either acoustic or con- 
vective) . The time delays to these negative peaks were removed and 
Fourier transforms of the translated correlations computed, to ob- 
tain that part of the phase shift not due to time delay. The re- 
sults indicated, in both cases, a 180° phase shift, independent of 
frequency. Similarly, the transfer functions between the com- 
bustor pressure and tailpipe pressures were measured, and their 
magnitudes were proportional to frequency squared. 

What we have seen, then, is that in the range of frequencies 
where the fluctuating combustor pressure is related to far-field 
noise, the tailpipe pressures behave approximately as the second 
time derivative of the combustor pressure. 

In a similar manner, we can now show that in this frequency 
range (0-240 Hz), pure propagation exists within the tailpipe. 

The broadband (0-1600 Hz) cross-correlation between the fluctuating 
pressures at the two tailpipe stations was shown in figure 12. The 
low pass filtered version of this correlation is shown in fig- 
ure 24. Here the even symmetry properties with respect to a posi- 
tive peak at a delay time corresponding to the acoustic delay pre- 
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The compressed version of this correlation was translated to 
remove this time delay and then Fourier transformed on the DSP . 

The phase part of the result, shown in figure 25, exhibits essen- 
tially a zero degree phase shift, independent of frequency. 

(The '"spike’’ is sixty cycle noise.) The amplitude of the trans- 
fer function between these two pressures, as measured directly by 
the DSP is shown in figure 26, for the low frequency range of in- 
terest. In contrast to the frequency squared relation which ex- 
ists between the combustor pressure and the tailpipe pressures, 
here the amplitude of the transfer function is flat and relatively 
independent of frequency. These two properties, 2 ero degree phase 
shift (after time delay removal) and flat transfer function, com- 
bined with a delay in the time domain corresponding to acoustic 
propagation, indicate pure acoustic propagation. 

We can now complete the picture by examining the phase and 
amplitude relationships between the combustor pressure and the 
far-field acoustic pressure, and between the tailpipe pressures 
and the far-field acoustic pressure. The corresponding . ow-pass 
cross— correlations, phase shifts with and without time delays, 
and transfer function amplitudes are shown in figures 27 
through 35. 

The combustor to far-field relations (for the 120° far-field 
angle) are shown in figures 27 through 29 and here we see a 180° 
phase shift (fig. 28(b)), frequency squared relationship 
(fig. 29). The far-field acoustic pressure, then, is also 


behaving as the second time derivative of the fluctuating pressure 
in the combustor. The upstream tailpipe to far-field relations and 
the downstream tailpipe to far-field relations are shown in figures 
30 to 32 and 33 to 35, respectively. In both cases, a zero degree 
phase shift and flat amplitude transfer function prevail, and these 
are indicative of pure time delay. 

The picture which emerges then, is the following: For the 

frequency range where the far-field acoustic pressure is coherent 
with the fluctuating combustor pressure, the far-field acoustic 
pressure, as well as the pressures within the tailpipe, behave as 
the second time derivative of the combustor pressure. This be- 
havior is consistent with a combustion noise generating mechanism 
wherein the fluctuating pressure due to entropy fluctuations acts 
as a source for acoustic energy, as discussed in Chapter 2. Within 
the tailpipe and between the tailpipe and far-field, for the same 
low frequency range, the evidence indicates pure acoustic propaga- 
tion. The nozzle, then, is acting as a transmission line for the 
acoustic energy generated in the combustor, at least for frequen- 
cies up to about 250 Hz. 

4.4 Variation of Data with Engine Operating Conditions 

and Far-Field Microphone Data 

The results presented so far in this chapter have been limited 
entirely to an engine operating speed of 43 percent of maximum, 
and a far-field microphone angle of 120°. These data, however, are 
typical of the results obtained at all microphone angles (40° 



through 160°) and at operating speeds from 30 percent (less than 
ground idle conditions) through 60 percent (approach conditions) of 
maximum fan speed. The conclusion, therefore, that there is a mea- 
surable contribution to overall engine noise from the combustor <>' 
holds true throughout the forward and rearward arcs of the engine, 
and for a substantial portion of the static operating conditions. 

The spectral range of this contribution, with only minor variations, 
peaks near 125 Hz and is limited entirely to frequencies below 
about 250 Hz. 

For the purpose of presenting supporting evidence for this, the 
results for the pair-wise relations between combustor pressure and 
far-field acoustic pressure at a forward arc angle of 60° are shown 
in figures 36 through 38. The coherence function, amplitude trans- 
fer function, and time delay removed phase shift are contained in 
figures 36, 37, and 38, respectively. The coherence, though not as 
strong as at the 120° microphone is seen to be limited to frequen- 
cies below 250 Hz, although the spectral peak is not as distinct. 
This coherence function is shown in appendix C to produce an un- 
usual directivity pattern for combustion noise at certain frequen- 
cies. Unfortunately, as is shown in appendix A, its reduced nu- 
merical value at some frequencies decreases the statistical 
confidence in the data and establishes a relative broad uncertainty 
band in the resulting coherence spectrum. Figure 37 shows a trans- 
fer function which, although not as clean as those shorn earlier, 
still exhibits a frequency squared behavior. Finally, the time 


delay removed phase shift, figure 38, indicates a phase shift of 
approximately 180°, over the frequency range of interest. Here 
again, however, the result is not as clean as at the 120° angle, 
but the trend is obvious. All these results at the 60° (forward 
arc) far-field angle, are for the same operating conditions as 
earlier in this chapter: 43 percent of maximum speed. 

Figures 39 to 42 show similar results between combustor and 
far-field at a different operating condition, 30 percent speed, but 
at the same rearward arc position of 120° shown earlier. The co- 
herence, amplitude transfer function, and zero time delay phase 
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shift are depicted in figures 39, 40, and 41, respectively . The 
coherence result shows essentially the same low frequency corre- 
lating region seen in the previous two sets of data. 

The amplitude transfer function (fig. 40) again shows a fre- 
quency squared behavior, and the zero-time delay phase shift 
(fig. 41) is near 180° at all frequencies within the range of in- 
terest. 

4.5 Discussion of the Diagnostic Results 

As indicated earlier, most previous investigators studying core 
noise have ~ound the dominant spectral peak to occur in the range of 
400 to 500 Hz. The present Investigation indicates a lack of any 
coherent relationship between combustor pressure and far-field 
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The tone at 380 Hz in fig. 39 was also observed in the com- 
bustor pressure spectrum and at almost all far-field microphones. 

It only occurred at 30 percent engine speed, which is an off-design 
operating point. 


acoustic pressure beyond about 250 Hz , Some question arises, there- 
fore, as to the validity of the coherence measurements. 

Reference 18, for example, presents the results of coherence 
measurements between the pressure measured in a gas turbine combustor 
installed in a specially designed test apparatus (i.e., non-engine 
installation) and far-field acoustic pressure. The authors distin- 
guish between acoustic pressures in the combustor and non- 
propagating hydrodynamic pressure fluctuations (i.e., pseudosound) . 

It is then pointed out that the presence of this pseudosound would 
serve to reduce the pressure coherence between the combustor and 
the acoustic far-field. The conclusion in reference 18 is that this 
reduced coherence Is not necessarily indicative of a correspondingly 
reduced combustion noise contribution to the far-field. Based on 
these results, therefore, it may be argued that this pseudosound 
"contamination” is responsible for the lack of any pressure co- 
herence between combustor and far-field beyond 250 Hz. 

It is the view of the present author, however, that in a source 
region this distinction between acoustic pressures and non-acoustic 
pressures loses, its meaning. In the present work the causal re- 
lationships between the fluctuating pressure in the source region 
and the acoustic far-field are being investigated. As such, we 
are interested in the entire pressure field in the combustor, what- 
ever its nature, and not just those pressures which can be specifi- 
cally identified as being locally acoustic in their character. The 
issue of turbulence "contamination" in the combustor pressure mea- 


surement, therefore, is not meaningful in the present context. 

The strongest supporting evidence for the spectral results re- 
ported in this chapter may be found in reference 33. Here the 
authors reported the results of a combustion noise investigation by 
far-field acoustic measurements alone conducted on a small aircraft 
auxiliary power unit (APU) . The APU is basically a small gas tur- 
bine engine used for auxiliary power generation while the aircraft 
is on the ground. As such, there is no fan or convergent nozzle to 
generate additional noise. 

The results reported in this reference, with the turbine re- 
moved, clearly indicate a combustion noise peak at 125 Hz. Further- 
more, at 250 Hz the measured sound power level is down approximately 
10 dB.^ This is consistent with the results to be shown in the next 
chapter, where the combustor-associated far-field power level spec- 
tra are presented. 

Finally, the authors of reference 33 indicate that this peak 
frequency remained fixed for all combustors tested in the small gas 
turbine class, even the reverse flow configurations such as inves- 
tigated here. It will be recalled from the description of the YF-102 
in Chapter 3, that the core of this engine is basically a small gas 

^The data in ref. 33 are reported in full octave bands, and the 
actual results show about a 7 dB reduction in combustor sound power 
between 125 and 250 Hz. The data reported herein is for constant 
bandwidths. The conversion from proportional bandwidth data to con- 
stant bandwidth data can easily be accomplished by subtracting 
10 dB/octave (^3 dB/octave) from the proportional bandwidth data. 

This conversion holds regardless of the precise width of the constant 
bandwidth data. 
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turbine in comparison with the cores of the turbofan engines in com- 
mercial service, and on which much of the combustion noise testing 
has been done. It would appear, then, that viewed with respect to 
the results reported in reference 33, the present results are rela- 
tively consistent. 

The results of reference 33 also prove useful in interpreting 
a portion of the diagnostic measurements which have not been pre- 
sented yet. In the previous sections of this chapter, all the co- 
herence results reported were with respect to the combustor pres- 
sure. Pressure coherence measurements were also made, however, 
between each of the tailpipe measuring stations, as well as between 
the tailpipe stations and the acoustic far-field. The results are 
shown in figures 42 to 44. Figure 42 shows the pressure coherence 
between the two tailpipe stations. Figures 43 and 44 show the co- 
herence between the upstream tailpipe station and the 120° far- 
field microphone, and between the downstream tailpipe station and 
the 120° far-field microphone. These data are for an engine speed 
of 43 percent of maximum, as before. 

These data show the same region of relatively high coherence 
in the 0 to 250 Hz region as prevailed with respect to the combustor. 
Now, however, there are additional broad regions of lower, but not 
insignificant, coherence at frequencies up to about 1 kHz. An at- 
tempt was made to analyze these higher frequencies regions with the 
same technique as was used for the low frequency combustor associ- 
ated noise. No simple result, however, emerged. 
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Figure 45(a), for example, shows the 250 to 1000 Hz band-pass 
filtered correlation between the downstream tailpipe pressure and 
the 120° far-field microphone. The correlation is seen to oscillate 
rather strongly at a delay time of about 87 msec, the acoustic 
travel time between the engine and the far-field microphone. Some 
of this oscillation, or "ringing" is characteristic of the correla- 
tion functions of band-pass random signals and gets worse as the 
pass band narrows , or the center frequency increases (26) . To some 
extent, then, it is an unavoidable result of the filtering process. 
The result, however, is that it is extremely difficult to determine 
which peak (positive or negative) , or which zero crossing corre- 
sponds to the delay time between the two signals. 

As described earlier, to obtain the zero time delay phase shift 
by Fourier transforming a translated cross-correlation function, it 
was necessary to recompute the correlation at a slower sampling rate. 
This is necessary to increase the percentage of useful points in the 
resulting frequency domain representation, and to decrease the reso- 
lution bandwidth to minimize bias error. This correlation compres- 
sion, of course, produces reduced time domain resolution. 

The result of recomputing this correlation at a slower sampling 
rate is shown in figure 45(b). Clearly, a substantial amount of de- 
tail has been lost. The DSP used in this investigation permitted 
continual visual monitoring of the phase roll on a CRT while this 
correlation was translated in time and Fourier transformed. No 
identifiable phase information was found through a complete rota- 




tiera of this correlation. In a similar manner, the transfer func- 
tion in this pass hand did not exhibit any simple properties which 
could be identified with either pure time delay (flat), surface in- 
teraction noise (proportional to frequency) , or combustor related 
noise (proportional to frequency squared) . 

Nevertheless, figures 43 and 44 demonstrate that a broad range 
of frequencies beyond 250 Hz is propagations! and is contributing 
to overall engine noise. Alternatively, these higher frequency re- 
gions of pressure coherence within the tailpipe and between tailpipe 
and far-field could result from external noise, especially jet noise, 
propagating back into the tailpipe. This may be ruled out in view 
of figure 12, which shows the correlation between the two tailpipe 
stations to have a single, positive peak, indicating downstream 
travel only. 

The results reported in reference 33 suggest a possible ex- 
planation. In addition to reporting data on the APU without the 
turbine, the authors also show data with the turbine installed. For 
this case, there is an additional region of core noise extending to 
beyond 1 kHz.^ The additional acoustic energy must, then have the 


■’As shown in ref. 33, this added noise with the turbine in place 
extends to about 10 kHz. The data is shown in full octave bands. 

The high frequency information, say beyond several kHz, is likely 
associated with the rotational speeds of the turbine. It does not 
show up as sharp peaks , or tones , because of the very wide bandwidth 
(approximately 2800 Hz at a center frequency of 4 kHz) of an octave 
analysis. In the region we are concerned with, less than 1 kHz, the 
data in ref. 33 are either flat or increasing slowly. This would 
correspond to a flat or slowly decreasing spectrum for the constant 
bandwidth analysis used here. 
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turbine as its source, and this appeal’s likely in the present case. 
Whether it is a result of surface interaction as the gases flow 
across the turbine surfaces, or the so-called "indirect" combustion 
noise (see ref. 34, for example) caused by entropy spots convecting 
through the turbine expansion passages is not evident. Neverthe- 
less, it may certainly be classified as core noise. However, the 
emphasis in the present study is on direct combustion noise, and it 
will not be further investigated. 

An additional point for discussion involves the possibility 
that at least some of the combustor-associated noise results from 
the interaction of the very turbulent combustion gases with the sur- 
faces of the combustor. The actual surface geometry of a combustor 
is rather complex with numerous holes for the admission of compressor 
discharge air, and the promotion of more efficient combustion. In many 
studies of gas turbine combustion noise, baseline tests are con- 
ducted in which heated or unheated air is passed through the com- 
bustor, but without the addition of fuel, so that no combustion ac- 
tually occurs. The mass flow rates used are normally the same as 
would exist with combustion, and these types of tests are generally 
performed in specially designed test apparatus (i.e., no engine). 

These tests are normally conducted to determine that part of 
combustion noise which is flow-related only; that is, scrubbing 
noise. As shown in reference 25, however, and briefly described in 
Chapter 2, if the flow passes over a finite surface (i.e., one with 
a leading or trailing edge) , the correlation between source pressure 
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and far-f ield acoustic pressure should have an anti-symmetric shape. 
It should be an odd function, and the associated phase and an^litude 
relations should be first derivative in nature: 90° phase shift and 

transfer function amplitude linearly proportional to frequency. The 
results shown here do not exhibit these properties. Consequently, 
it must be concluded that any scrubbing noise when actual combustion 
is present is negligible compared to the noise directly associated 
with combustion. 

The final remarks to be made in this chapter address the issue 
of a very important requirement for proper interpretation of cor- 
relation and coherence measurements, and the meaningfulness of a 
transfer function: the system being studied with these random data 

analysis techniques must be linear. This, of course, has been an 
implicit assumption in all that has preceded. If the acoustic 
propagation mechanism between any pair of points being investigated 
is non-linear, then the interpretation of the measured transfer func- 
tions is, at very least, questionable. Furthermore, non-linear 
propagation would certainly serve to reduce the value of the co- 
herence between pressure measurements. This may be advanced as a 
possible reason for the virtual vanishing of any pressure coherence 
between the combustor and far-f ield acoustics beyond 250 Hz. 

It is likely that if a nonlinear propagating mechanism were 
present, it would result from high pressure levels present in the 
combustor. The spectrum of the fluctuating combustor pressure up 
to 1 kHz, at 43 percent operating speed, is shown in figure 46. The 
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level at 250 Hz, where the coherence has essentially gone to 2 ero 
(see fig. 16) is about 118 dB (approximately 0.0023 psi, rms) . This 
is certainly not excessively high, even by acoustic standards. It 
does not seem reasonable that this pressure level is sufficiently 
high to result in such predominantly non-linear propagation as to 
decorrelate the pressure signals . Nevertheless , without a detailed 
knowledge of the precise propagation mechanism the presence of non- 
linear effects must be admitted as at least being possible. In view 
of the previously cited evidence from reference 33, however, the 
assumption of negligible nonlinear effects is still felt to be valid. 




5. MEASUREMENT OF COMBPSTION NOISE 


5.1 Introduction 

The previous chapter described the results of a series of diag- 
nostic pressure measurements conducted on an AVCO-Lycoming YF-102 
turbofan engine. The measured ordinary coherence functions between 
the fluctuating pressures at several locations within the core of the 
YF-102 and between these pressures and the acoustic pressures in the 
far-field were used only as a frequency domain guide for measuring 
filtered cross-correlation functions . The shapes of these filtered 
correlation functions were then used to make a qualitative judgement 
that the local fluctuating combustor pressures were acting as 
sources for low frequency acoustic waves propagating through the 
tailpipe and out to the far-field. The phase and amplitude informa- 
tion obtained from the appropriate Fourier transforms of these cor- 
relations were shown to be consistent with these conclusions. In 
this chapter, the ordinary coherence function will be used to obtain 
quantitative results on the far-field acoustic contribution of com- 
bustor generated noise. 

At any given frequency, for the idealized case of a perfectly 
linear system with one noise-free input and one noise-free output, 
the ordinary coherence function between input and output has a value 
of unity. Conversely, if the input and the output of the system are 
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completely unrelated, then the ordinary coherence function must have 
a value of zero. The ordinary coherence function may have a value 
between zero and one if one or more of the following three cases 
exist (26) : 

(a) Contaminating (i.e., noncorrelating) "noise" is present in 
the measurements of input or output. 

Cb) The output is due to inputs other than, or in addition to, 
the input being measured. 

(c) The system relating the output to the input is nonlinear. 

In the present context, the fluctuating pressure in the combus- 
tor is considered to be the "input" and the far-field acoustic pres- 
sure is considered as the "output." The "system" is the generating 
mechanism and propagation path between the combustor and the far- 
field microphones , including any path within the engine through which 
the pressure signal travels. It has been assumed throughout the 
present investigation that this system is, in: fact, linear. A brief 
argument supporting this view was presented in the previous chapter. 
It was also argued that the "pseudosound" pressures within the com- 
bustor could not be considered as a signal contaminant when investi- 
gating casual relations with the combustor as a source region. On 
this basis, then, the value of the ordinary coherence function may 
be considered a measure of the far field acoustic energy which is due 
to the fluctuating pressure within the combustor. Conversely, the 
quantity one minus the coherence function is a measure of the acous- 
tic energy not due to the fluctuating combustor pressure. It must, 
therefore, be a measure of the far-field noise due to "inputs" such 
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as fan noise, jet noise, and other internal noise sources which do 
not correlate with the fluctuating combustor pressure. In equation 
form, then: 

C(f) = Yc F ( f ) G FF ( £ ) f “ frequency (30) 

Here, C(f) is the combustor coherence spectrum; that is, it is that 

portion of the far-field spectrum which is combustor related, and is 

therefore the combustion noise contribution to the far-field. 

Y'L(f) is the coherence function between fluctuating combustor pres- 

sure and far-field pressure; and G (f) is the far-field acoustic 

r c 

spectrum. 

We note, however, that it is possible that the single point 
pressure measurement within the combustor is not adequate to charac- 
terize the entire combustor source region. That is, there may be 
several (or perhaps, many) independent source regions within the 
combustor, each contributing to the far field acoustics. If such is 
the case, then the combustor coherence spectrum as computed from 
equation (30) will not include the contributions from the other com- 
bustor source regions and hence will be an underestimate of the true 
combustion noise contribution. (Note that this is an extension of 
case (b) above.) Evidence will be presented shortly, though, that 
this single point combustor pressure measurement is sufficient to 
make at least a reasonably accurate estimate of the actual combus- 
tion noise contribution to the far-field, and is entirely sufficient 
to characterize the spectral content and directivity of the combustor 
associated far-field noise. 
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5.2 Statistical Considerations 

Since the techniques used here to measure the combustion noise 
utilize the tools of random data analysis, statistical errors, which 
are an inevitable byproduct of such techniques* are an important con- 
sideration. A discussion and computation of the confidence limits 
associated with the results to be shown are presented in appendix A. 
However, one particular source of error in the measurement of the 
coherence functions is related to the question of whether or not the 
single point combustor pressure measurement is sufficient to char- 
acterise the entire combustor source region, and so will be dis- 
cussed here. 

As shown in reference 30, and briefly mentioned in Chapter 3, 
an important bias error occurs when examining pairwise relations be- 
tween signals with significant time delay between them. The proper 
way to eliminate, or at least minimize, this bias error is to remove 
this delay before processing, by shifting the appropriate signal in 
time. This shifting can be readily accomplished if the digitized 
data is being processed on a user-programmable digital computer. 
However, the digital signal processor used in this study did not 
have a provision for implementing this precomputation delay when 
computing coherence functions. Consequently this source of bias 
error could not be removed from the raw measurements. 

In reference 30, the bias error due to the time delay is shown 
to be given by 
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Y 2 (f) = Y 2 (f)(l - (31) 

where 

Y 2 (f) is the coherence function estimate (i.e., the measured value) 

2 

Y (f) is the true value of the coherence function 
t is the delay time between the two signals 

and 

T is the processor memory period (i.e., the sampling period 

for a single record) (See appendix A.) 

In appendix A it is shown that for the computation scheme uti- 
lized by the DSP, the memory period is inversely proportional to the 
analysis frequency range f^, and is given by T = 500/f^. With 
this brief discussion* we can now return to the issue of the adequacy 
of a single point pressure measurement within the combustor in com- 
puting the combustor coherence spectra. 

5.3 The Combustor as a Single Source Region 

The first point to be made here is that one relative straight— 
forward way of determining whether or not the combustor consists of 
a multitude of individual sources is to make a series of space-time 
cross-correlation measurements around the periphery of the combustor. 

The zero time delay value of these correlation functions can then be 
cross-plotted as a function of circumferential separation with re- 
spect to some fixed position and the result will be a two point cor- 
relation with reparation distance as the independent variable, rather 
than time. From this result we could then estimate the circumferen- 
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tial correlation length of a single source region, and hence obtain 
the number of independent source regions within the combustor. How- 
ever, only two fixed pressure measuring stations were available 
within the combustor and this approach could not be used. 

Alternatively we can examine the pressure coherence between the 
two measuring stations available. If this results in a uniformly 
high coherence over the frequency range of interest (approximately 
0 to 250 Hz) , then this would at least suggest that a 90° sector is 
the minimum size of the correlating region and that, at most, there 
are about four independent source regions. This coherence function 
between the pressure signals 90° apart in the combustor is shown in 
figure 47. The result shows that uniform coherence does not prevail. 
The coherence between the two is relatively uniform and high frcm 
about 0 to 100 -Hz or so, but drops substantially beyond 100 Hz. This 
result is consistent with the notion that low frequencies tend to 
correlate over longer distances than high frequencies. It suggests 
that up to about 1(00 Hz, the correlating (i.e., source) region ex- 
tends to at least 25 percent of the combustor annulus. At frequen- 
cies beyond about 100 Hz, however, the results seem to indicate that 
the correlating regions are smaller and therefore more numerous. 

This issue can be at least partially resolved by reexamining 
the coherence function between the pressure in the combustor and the 
acoustic pressure in the far-field. This function was originally 
shown in figure 16. There, however, we were concerned with diag- 
nostic measurements only and the numerical values of the coherence 
were not of immediate concern. To illustrate more graphically the 
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fact that there is virtually no acoustic contribution from the com- 
bustor to the far field beyond about 250 Hz, a relatively broad anal- 
ysis range of 1500 Hz was chosen to present the data. The memory 

period for the DSP at 1500 Hz is T = 500/f = 500/1500 - 0.333 

m 

seconds. From figure 27(a), the time delay, t, between combustor 
and far field is about 87 msec. So, from equation (31), the bias 
error due to time delay is: 



The coherence estimate in figure 16 then is only about 54.5 percent 
of the true value. 

To show the effect of this time delay bias, the coherence be- 
tween combustor and 120° far-field pressures was recomputed at a 
lower frequency range (at the expense of increased random error, as 
is shown in appendix A) . The result, for an f of 400 Hz is shown 
in figure 48. Now the peak coherence has risen to a value of about 
0.62, in contrast to the value of 0.33 in figure 16. We could re- 
compute this coherence function once again at a still lower fre- 
quency of 200 to 250 Hz, the maximum frequency of interest. Here, 
however, the random error gets too large (see appendix A), so a cal- 
culated bias correction will be applied to the results of figure 48. 
For this case, the memory period of the DSP is T = 500/400 = 1.25 
seconds. (The time delay, t, is still 87 msec.) 
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Y 2 (f) 



0.62 

0.869 


0.72 


The true value of the coherence, then, is nearer to 0.72. So, 
at the peak frequency, approximately 72 percent (i.e., within 
1.5 dB) of the far-field signal can be accounted for by a single 
point pressure measurement in this combustor. Considering the fact 
that at least part of the far-field noise at this frequency must 
contain contributions from the other engine noise sources, the as- 
sumption that a single point pressure measurement in the combustor 
is sufficient to characterize the source region appears reasonable. 
This argument applies only to frequencies near the peak. However, 
again considering the notion that low frequency information should 
correlate over larger distances than high frequency information, 
the assumption is probably valid at frequencies below the peak also. 


when even, larger correlation volumes could be expected. 

We are now left only with the question of the reduced coherence 
at frequencies between about 125 Hz and 200 to 250 Hz. With the in- 
formation obtained to this point, no firm conclusions can be reached 
as to whether this results from multiple, independent source regions 
within the combustor at the higher frequencies (i.e., above 125 Hz) 


or is a true indication of the combustor-associated far-field noise. 


The overall combustor noise power levels in the far-field, which 
are presented in a later section, will be seen to be somewhat low 
when compared with most of the empirical predictions, although not 
excessively low. The possibility is admitted, therefore, that 


additional measuring stations within the combustor may be required 
to account for all the source regions within the combustor. Their 
exact number, however, could not be determined without a comprehen- 
sive analysis or experimental program to completely characterize 
the combustor as a source region. 

The possibility of different azimuthal regions within the com- 
bustor contributing to different spectral ranges in the far-field 
acoustics may be ruled out on simple geometric symmetry considera- 
tions. Figure 49, which shows the coherence function between the 
second combustor probe and the 120° far-fleld microphone is adequate 
evidence for this conclusion. The spectral range of coherence, and 
the shape of the coherence function, is virtually identical to that 
of the in-line combustor probe, shown in figure 16. 

5.4 Combustor Coherence Spectra; Radiation Patterns 

From equation (30), the combustor coherence spectrum is defined, 
for the present case, as: 

c(f) - Yj F (f)G EF (£) 

This can be computed directly through a frequency by frequency 
multiplication of the measured pressure coherence function between 
combustor and far-field and the measured auto-spectrum of the far- 
field acoustic pressure^". This multiplication could be implemented 

■^The result of this multiplication is often referred to as the 
coherent output power in literature on random data analysis. How- 
ever, to avoid confusion with the physical quantity of acoustic 
power, which will be calculated later, the term combustor coherence 
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directly on the DSP. 

The data to be shown were all computed with an analysis range, 
f^, of 400 Hz. In addition to the bias error in the estimate of y 
due to time delay, as described earlier, there are random and bias 
errors in the estimate due to finite sample records. As shown in 
appendix A, for the computational scheme utilized by the DSP, the 
criteria for minimizing all these errors are in conflict . The 
choice of 400 Hz, therefore, represents a compromise. 

The bias error due to time delay for f - 400 Hz was calcu- 
lated earlier, and found to underestimate the true coherence by a 
factor of 1/0.869. Expressed in decibels this error is seen to be 
4 = 10 log 0.869 - -0.65 dB. The computed combustor coherence 
spectra, then, will be biased low by about 0.6 dB. This small cor- 
rection is not made in the data to be shown in this section, but is 
accounted for when the acoustic power results are presented. Con- 
fidence limits for the combustor coherence spectra are described in 
appendix A and are tabulated in appendix B. 

Additionally before presenting the data, we note that correc- 
tions to the far-field spectra due to atmospheric absorption have 
not been made. At the low frequencies of interest here, this cor- 
rection is entirely negligible (35) . The data have also not been 
corrected to free-field to account for the pressure doubling due to 
the ground mounted microphones (+6 dB) . Again, however, these 


spectrum will be used here. This is in keeping with the idea that 
it represents that portion of the far-field acoustic spectrum which 
is coherent with the fluctuating pressure in the combustor. 
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corrections were applied to the acoustic power results, shown later. 
Also, no corrections were made to standard-day conditions. The data 
to be presented in this section, therefore, are "as-measured . " 
Finally, all the data to be shown in this section have been normal- 
ized to a spectral density basis, so that the results are presented 
in dB/Hs. 

Some representative far-field sound pressure level spectral 
densities at angles of 60°, 120°, and 160° are shown in figures 
50(a), (b), and (c) (dashed curves), respectively, for an engine 
operating speed of 43 percent. The resulting combustor coherence 
spectra at these angles are shown by the solid curves on the same 
figure. These were obtained by a direct frequency by frequency mul- 
tiplication of the far-field spectra by the appropriate coherence 
functions. The primary observation to be made is that because the 
far-field acoustic spectra are relatively flat over the relevant 
frequency range, the combustor coherence spectra peak approximately 
at the same frequency as the coherence functions themselves, around 
125 Hz as shown earlier. Also, as first suggested from the co- 
herence functions, the results indicate that there is virtually no 
combustor associated noise at frequencies above about 200 to 250 Hz. 

Results similar to above were obtained for each of the far- 
field microphones (40° through 160°) and for each of the engine 
operating conditions through 60 percent. Since the resultant com- 
bustor coherence spectra were all normalized on a spectral density 
basis, the area under any given combustor coherence Bpectrum repre- 
sents the overall sound pressure level (OASPL) , over the frequency 
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range 0 to 400 Hz (although the contribution beyond about 200 to 
250 Hz is negligible, as shown above.) . However, because the DSP was 
not interfaced with a computer this integration was performed 
manually for each combustor coherence spectrum. Eleven points on 
each spectrum, in intervals of 20 Hz, from 40 to 240 Hz were selected 
for the summation. These eleven points were anti-logarithmically 
summed. The result was then multiplied by 20 Hz (i.e., 

10 log 20 = 13 dB was added) to account for the effective 20 Hz band- 
width. The final result for each combustor coherence spectrum, then, 
represents an approximation to the OASPL up to a frequency of 240 Hz. 
For smoothly varying spectra, the approximation is adequate. 

This procedure was repeated for each microphone angle. The 
result may be plotted as a function of angle to produce the radia- 
tion pattern of the combustor associated far-field noise. The re- 
sults are shown in figure 51, for three engine power settings of 
30, 43, and 60 percent. As would be expected, there is a slow in- 
crease in combustor associated far-field noise at all angles as en- 
gine power is increased. The radiation patterns, however, remain 
constant, with a slight peak near 120°. The peak angle is con- 
sistent with that found by other investigators (ref. 35, for ex- 
ample) . Similar results were found to prevail at the intermediate 
engine power settings of 37 and 50 percent, but are not shown on 
figure 51 in the interest of clarity. 

Because the radiation patterns remain invariant with engine 
power setting, they may be normalized, for convenience purposes, 
into a single directivity index curve. The directivity index is 
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defined as the difference (In decibels) between the measured sound 
pressure level at a point and the sound pressure level that would 
exist at the same point from a simple monopole source emitting the 
same acoustic power. The simple-source sound pressure levels were 
computed by area weighting the combustor-associated OASPL, as de- 
scribed in reference 36. The computations were performed by using 
a slightly modified version of a standard computer program in use at 
the Lewis Research Center for handling acoustic data, also described 
in reference 36. The result, arithmetically averaged over the five 
engine operating speeds between 30 and 60 percent, is shown in fig- 
ure 52 (symbols) . 

For comparison purposes, the directivity index as obtained from 
the empirical prediction schemes outlined in references 19, 32, 
and 37 are also shown. The present results are seen to be in re- 
markably close agreement with the empirical prediction scheme of 
(37). The recommended prediction procedure of reference 19 appears 
to overpredict the results found here at forward and rearward angles 
near the engine axis. The procedure of (32) produces results which 
somewhat overpredict the measured data at the rearward angles. 
However, in the range between about 60° and 130°, the results found 
here agree reasonably well with all three prediction schemes. 

5.5 Combustor Associated Far-Field Acoustic Power 

The combustor associated far-field acoustic power spectra and 
total power may also be obtained from the combustor coherence spec- 
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tra. For each of the 11 frequencies selected from the combustor co- 
herence spectra, a point on the combustor far-field acoustic power 
spectrum was computed for a given engine power setting. The basic 
equation used was (36) : 



Here, p and c are the ambient density and speed of sound, respec- 
2 

tively; p F . is the mean-squared pressure at the selected frequency 
f, i is the summation index corresponding to each microphone loca- 
tion and the AA^ are finite incremental areas on which the corre- 
2 

sponding p f * are presumed constant* Any set of consistent units 
r ? i 

may be used for p, c, and p. Again, a standard computer program 
described in reference 36, was adapted for implementing the com- 
putation. 

The results, on a spectral density basis, (dB/Hz, referenced to 
-13 

10 W) are shown in figure 53, for three engine speeds of 30, 43, 
and 60 percent of maximum. The spectral shapes are seen to vary 
somewhat as operating speed increases with the primary differences 
occurring at frequencies beyond the peak. There appears to be a 
trend towards a slightly increasing peak frequency with engine 
operating speed, but this may be an artifact of the relatively large 
bandwidth (20 Hz) chosen for the numerical integration. 

-13 

The integrated overall sound power level (OAPWL) (re 10 W) 
are shown adjacent to each curve in figure 53 (corrected for time 
delay bias) . Table I presents a comparison of the present results 
with the empirical prediction schemes of references 19, 32, and 37. 
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The coherence results are seen to agree reasonably well with the pre- 
diction of reference 19, but are considerably less than those of 
references 32 and 37. This suggests the possibility of the inade- 
quacy of a single point combustor measurement in characterizing the 
source region, as discussed in section 5.2. 

Table I 


Overall Power Level 
(dB, re: 10 -13 W) 


Engine 
Speed , 
percent 
of max. 

Ref. 19 

Ref. 32 

Ref. 37 

Coherence 

results 

30 

120.9 

127.0 

128.7 

120.0 

43 

122.7 

131.5 

133.2 

123.5 

60 

126.3 

136.0 

137.7 

126.2 


The final comparison to be made is the peak frequency found by 
the coherence results here with the peaks predicted, for this engine 
by the procedures outlined in references 19, 32, and 37. For the 
most part these are independent of operating condition (as found 
here) within about 5 Hz. The results are: 


Ref. 

Predicted 
peak fre- 


quency , 


Hz 

19 

260 

32 

400 

37 

330 

Present 

125 


At engine speeds beyond 60 percent of maximum the numerical 
value of the coherence function between the combustor pressure and 
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the far-field pressure diminished significantly. This is shown in 
figure 54 for an engine speed of 75 percent, the next highest power 
setting which was tested. The low numerical value of the coherence 
results primarily from the rapid increase in jet mixing noise which, 
as will be recalled from the INTRODUCTION, increases as the eighth 
power of jet velocity. Beyond 75 percent speed, no measurable co- 
herence between the pressure in the combustor and the far-field was 
observed. This does not mean, of course, that the combustor is no 
longer contributing to overall engine noise, but that its propor- 
tional contribution at these higher engine speeds is extremely small. 

In principal this contribution, however small, can be extracted 
through the coherence function. However, as shown in appendix A, as 
the estimated value of the coherence function becomes small, the 
variance of the estimate becomes very large, for a fixed number of 
sample records. To reduce this variance requires the availability 
of large amounts of data* The 120 second record lengths used here 
are not sufficient to provide adequate statistical confidence when 
the numerical value of the estimated coherence function falls below 
about 0.05o This corresponds to a contribution of the combustor 
associated far-field noise some 13 dB below the overall level. As 
is shown in the tabulated confidence limits in appendix B, for a 
90 percent confidence limit of about ±3 dB, this would require ap- 
proximately five times more data than was used in this investigation, 
or record lengths of approximately 10 minutes. 
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Before concluding this chapter, two additional points regarding 
the numerical results found here should be discussed. Fi'rst, it is 
clear that there is a significant discrepancy between thf combustion 
noise spectral peak found here compared with that predic ;ed by the 
empirical procedures of other investigators. It should be pointed 
out that the data used to develop those empirical equations were ob- 
tained from engines and combustors quite different than the YF-102, 
with its reverse flow combustor (although the YF-102 operating and 
geometric parameters were used to arrive at the predicted peak fre- 
quencies shown earlier). Since it is not unreasonable to expect 
geometry to play a significant role in determining spectral charac- 
teristics, it is not surprising to find such discrepancies between 
the measured and predicted spectral peaks. Clearly, the investiga- 
tion of scaling relationships may prove a fruitful avenue for future 
research . 

Finally, we address the issue of the basic numerical accuracy of 
the data shown in this chapter. As pointed out in Chapter 3 the mea- 
surement errors, being associated entirely with the microphones, may 
be considered quite negligible when considering the accuracy of pre- 
sent day microphone and amplifier systems (sfel dB) . The largest area 
for inaccuracy results from the statistical errors associated with 
the data processing. A thorough discussion of these is presented 
in Appendix A. Briefly, it is shown that tue error in the integrated 
results (i.e., OASPL and OAPWL) are dominated by the errors near the 
corresponding peak spectral values, where the measured coherence is 
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quite high, resulting in relatively small errors. For example, at 43% 
of maximum speed, the 90% confidence limits on the calculated OAPWL is 
-1.7 dB to 1.0 dB, with virtually identical results at the other 
speeds. By contrast, the errors at the extremes of the spectrum are 
quite large due to low coherence. For 43% of maximum speed, the 90% 
confidence limit on the power level is about -3 dB to +2 dB at 40 Hz 
and about -10 dB to +4 dB at 240 Hz. Similar results prevail at the 
other engine speeds. 
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6. SUMMARY AMD CONCLUDING REMARKS 


6 . 1 Summary 

Chapter 1 . In chapter 1, the motivation for the experimental 
program described in this thesis was established. Specifically, it 
was indicated that recent evidence suggests that the reduction of 
fan and jet noise under flight conditions results in core associated 
noise becoming a significant contributor to overall turbofan engine 
noise. Furthermore, current research has shown the duct-burning 
turbofan to be a promising candidate for possible supersonic trans- 
port applications. In both cases above, the identification and 
measurement of combustion-associated far-field noise is required to 
characterize the overall engine noise. Such a step is necessary 
before further major reductions in aircraft propulsion system noise 
can be accomplished. 

Chapter 2 . In chapter 2, a review of Ligh thill’s governing 
equation for aerodynamic noise generation was presented. It was 
shown that in a combustion region, a slight modification in the 
basic source term of Lighthill’s equation enabled the fluctuating 
pressure in the combustor to be related to the far-field acoustic 
pressure. It was further shown that this relationship was uniquely 
determined by the phase and amplitude information contained in the 
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transfer function between combustor pressure and far-field pressure. 
Specifically, for the model proposed, the far-field acoustic pressure 
was related to the combustor pressure by a transfer function propor- 
tional to frequency squared, and a 180° phase shift. Furthermore, 
for the model proposed, it was shown that the corresponding time 
domain relation between combustor pressure and far-field acoustic 
pressure, the cross-correlation function, could be expected to have 
an evenly symmetric shape with a negative peak at a time delay cor- 
responding to the acoustic propagation time between the combustor 
source region and the acoustic observation point. 

These relations were contrasted to the corresponding information 
which would be expected for "scrubbing" noise (i.e., surface inter- 
action noise), or pure time delay. 

Chapter 3 . This chapter described the experimental hardware 
and instrumentation. Briefly, the engine on which the measurement 
program was conducted was an AVCO-Lycoming YF-102 turbofan engine 
with a bypass ratio of six to one. The fluctuating pressure mea- 
surements within the engine core were accomplished with specially 
constructed "semi-infinite" waveguide probes, with conventional 
microphones being used as the pressure transducers. The data ob- 
tained from these probes, as well as from an array of far-field mi- 
crophones, were off-line processed on a Fast Fourier Transform 
(FFT) digital signal processor (DSP) . The DSP was used to obtain 
the time and frequency domain relations described in chapter 2, as 
well as the ordinary coherence functions, between the fluctuating 
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pressures within the engine core and between these pressures and the 

far-field acoustic pressures. 

Chapter 4 . The results of a series of diagnostic measurements 
were presented in this chapter. The measured coherence function be- 

- 

tween the fluctuating combustor pressure and the far-field acoustic 

pressure indicated that the fluctuating combustor pressure was not ^ 

related to the far-field acoustic pressure at frequencies beyond 
about 25G Hz. Using this as a guide, low-pass filtered cross- 
correlation functions were computed between combustor pressure and 
tailpipe pressure, between combustor pressure and far-field acous- 
tic pressure, and between tailpipe pressure and far-field acoustic 
pressure. 

The shapes of these functions as well as the corresponding fre- 
quency domain information (transfer function magnitude and phase) 
obtained by Fourier transforming them after removing the time de- 
lays, were consistent with the behaviour modeled in chapter 2. 

Specifically, the transfer functions between combustor pressure and 
nozzle pressure, and between combustor pressure and far-field acous- 
tic pressure were found to have a magnitude approximately propor- 
tional to frequency squared. Similarly, the corresponding phase 
shifts were found to be approximately 180° independent of frequency, 
after time delay removal. By contrast, the transfer functions be- 
tween nozzle pressure, and between nozzle pressures and far-field 
pressure were found to have a magnitude independent of frequency 
(i.e., flat) and a phase of 0° after time delay removal. From these 
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results it was concluded the combustor was a source region for low 
frequency noise (less than 250 Hz) which propagated through the noz- 
zle and out to the far-field. 

Chapter 5 - Having identified the combustor as a source region 
for far-field noise in the previous chapter, its quantitative con- 
tribution was determined in chapter 5, The tool used for this was 
the ordinary coherence function between combustor pressure and far- 
field pressure. A frequency by frequency multiplication of this 
ordinary coherence function by the corresponding far-field spectrum 
produced the combustor coherence spectrum. Under the appropriate 
assumptions, the combustor coherence spectrum may be interpreted as 
the spectral contribution of the combustor to overall engine far- 
field noise. These spectra were obtained for each microphone angle 
and over a broad range of engine operating conditions. These re- 
sults were, in turn, used to obtain the directivity patterns of the 
combustor associated far-field noise, acoustic power spectra, and 
total acoustic power. 

The directivity patterns were found to be relatively independent 
of engine operating speed, and agreed well over most angles with two 
existing prediction schemes, and remarkably well over all angles with 
a third. The power spectra were found to peak near 125 Hz, rela- 
tively independent of engine speed, with greater amounts of energy 
at frequencies beyond the peak as engine speed increased. In all 
cases, however, there was virtually no combustor associated noise 
at frequencies above about 250 Hz. The total acoustic power was 
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found to agree very closely with that predicted by one of three 
existing prediction schemes with which it was compared, but some- 
what low when compared with the other two. The peak frequency of 
125 Hz was found to be approximately one to two octaves below the 
results obtained from all three prediction schemes. Reference evi- 
dence was cited, however, to support this result for gas turbine 
combustion noise peak frequencies for the size class of this engine. 

6 . 2 Concluding Remarks 

With respect to the specific problem of combustion noise from 
turbofan engines, clearly there is much work still required to en- 
able the development of techniques for its suppression to be ac- 
complished on a rational basis. For example, the design of liners 
or suppressors for combustion noise reduction obviously requires 
knowledge of its spectral content. As was shown in the previous 
chapter, however, the existing empirical schemes for predicting 
turbofan combustion noise characteristics are certainly not adequate 
for predicting spectral content. Yet, extensive parametric testing 
of suppressor hardware or new combustor designs in full scale oper- 
ational turbofan engines using the techniques described in this in- 
vestigation would be prohibitively expensive. Clearly, it would be 
desirable to implement such testing on less expensive, more adapt- 
able apparatus specifically designed for such purposes. However, as 
pointed out earlier, it is not yet clear how the results of combus- 
tion no.' se tests in such apparatus can be related or applied to 
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actual engines. Here, it would seem, the techniques and procedures 
described in earlier chapters may be used or extended to enable a 
fuller understanding of combustion noise phenomena in such combustor 
rigs. 

Finally, as stated in the INTRODUCTION, one of the primary 
objectives of the present investigation has been to specifically 
identify combustion noise as a component of core associated far- 
fieid noise from a turbofan engine. As described in chapter 4, 
this identification has been made. Furthermore, in contrast to 
procedures which rely solely on far-field acoustic measurements, 
the identification of the combustor as a source for low frequency 
core associated noise is unambiguous. Of equal significance, how- 
ever, has been the demonstration that the techniques of random 
data analysis can be successfully applied to full scale, operational 
turbofan engines for both diagnostic and measurement purposes. This 
success, clearly, depends to a large extent on having a physical 
model of the acoustic generating mechanism available. The measure- 
ment techniques described herein, as well as the physical models for 
aerodynamic noise generation have been available for many years. It 
Is the bringing together or both in a unified experimental program 
with practical engineering objectives which the present author feels 
is the major contribution of this work. 
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Figure L - Turbofan Engine Noise Sources. 
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Figure Z - Source-field Geometry, 
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Figure 6. - Core Pressure Probe Locations. 
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Figure 8* - Schematic Of Core Probe Installation. 
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Figure 9. - Frequency Response For A Typical Acoustic Waveguide Probe. 
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Figure la - Yh -102 Turbofan Engine With Mounted Probes. 


Figure 11. - Test Arena Geometry. 
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Figure 15* - Coherence Between Combustor and Downstream Tailpipe Pressures (43% Speed). 



Figure 16* - Coherence Between Combustor Pressure and 120° Far-Field 
Acoustic Pressure (43% Speed)* 
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Figure 19. - Compressed Cross -correlation Between Combustor And Upstream Taitpfpe 
Pressures {42*. Speed), 



Figure 20. - Phase Shift Between Combustor And Upstream Tailpipe Pressures. Time Delay Removed t43-* Speed), 



Figure 21, - Transfer Function Magnitude Between Combustor And Upstream Tailpipe 
Pressures (43ft Speed). 


r.WSffvjfaA-. 



PHASE SHIFT, deg 


104 





PHASE SHIFT, DEG 





DELAY TIME, MSEC 

Figure 27* - Low-pass Filtered Cross-correlation Between Combustor Pressure And 12CP Far-fTeld Acoustic 
Pressure 143% Speed)* 
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Figure 28. - Phase Shift Between Combustor Pressure and 120° Far-Field 
Acoustic Pressure (43% Speed), 
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Figure 29. - Transfer Function Magnitude Between Combustor Pressure And 12CP Far 


field Acoustic Pressure (43% Speed), 
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Figure 3a - Low Pass Filtered Cross-correlation Between Upstream Tailpipe Pressure and 12QP Far-field Acoustic 
Pressure (43% Speed). 
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[gure 41. - Phase Shift Between Combustor Pressure and 120° Far-fleld Acoustic Pressure; Time Delay Removed 
(30% Speed). 
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Figure 42, - Coherence Between Tailpipe Pressures (43% Speed). 
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Figure 43. - Coherence Between Upstream Tailpipe Pressure And 120° Far- 
field Acoustic Pressure (43% Speed). 
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Figure 46, - Combustor Pressure Level Spectrum, 43ft Speed 



Figure 47, - Coherence Between Combustor Pressures; 
90° Azimuthal Separation (43% Speed), 
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Figure 48, - Coherence Between Combustor Pressure And 120° Far-field Acoustic Pres 
sure (43% Speed), 
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APPENDIX A 


STATISTICAL CONSIDERATIONS 


The purpose of this appendix is to discuss some of the statis- 
tical uncertainties in the data presented in this study. These 
statistical errors are an inevitable byproduct of random data anal- 
ysis techniques and, for the most part, cannot be entirely elimi- 
nated from the resulting data. 


Computational Scheme 


The digital signal processor used to obtain the results pre- 
sented in this s tudy^ digitized each sample record of data into 
1024 words at a sampling rate 2.048 times the highest frequency 
(f^) selected for analysis : 


f = 2.048 f 

sample m 


The processor memory period, or sample record length, therefore, 


is: 


_ Number of words/record 
- Number of words/sec 

= 1024/(2.048 f ) = 500/f , sec/record (Al) 

m m 

Since the number of words/record is fixed, the resolution bandwidth 
is 

B = 1/T = f /500 
e m 
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120 


The number of statistical degrees of freedom per sample record, 

therefore is (ref. 26) 

n - 2B T - 2(f /500) (500/f ) - 2 
e m m 


(A2) 


So, if N disjoint (i.e., independent) sample records are averaged, 
the total number of statistical degrees of freedom in the average is 


n = 2N 
tot 


(A3) 


For a total tape record length of 120 seconds, the maximum 
number of independent sample records is 

N = 120/T (A4) 


Bias Error Due to Time Delay 

In the usual analysis of signal pairs, the similarities, or 
mutual properties between the two signals do not coincide with one 
another until some time period has elapsed. In the present inves- 
tigation, of course, this is the acoustic propagation time between 
the combustor and the far-field, which was shown in chapter 4 to be 
about 87 msec. 

As pointed out in reference 30, failure to account for this 
time delay leads to a bias error in the coherence function, which 
becomes underestimated. The best procedure to eliminate or minimize 
this bias error is to delay the signal which occurs earlier in time, 
thus allowing the later signal to "catch up" before implementing any 
computation. The processor used to obtain the results shown in this 
paper, however, did not have such a precomputation delay feature 
when computing coherence functions. The bias error, therefore, was 


calculated and the correction applied to the overall power levels. 
From (30), the bias error due to time delay is given by 


vL(f) V / 


(A5) 


where 

is the estimated (i.e., measured) value of the coherence 
function between combustor and far-field, at the fre- 
quency, f 

2 

Y C p(f) is the true value of the coherence function 

t is the delay time between the two signals 

T is the processor memory period 

It is obvious from (A5) that to minimize this bias error, T 
should be made as large as possible. From equations (A3) and (A4) , 
however, to increase the number of statistical degrees of freedom, 
and therefore minimize random error (see below) , T should be made 
as small as possible. Any selection of T, therefore, in computing 
the coherence function, represents a compromise. For the results 
shown in chapter 5, f was chosen at 4Q0 Hz. So, from (Al) , 

T = 500/400 Hz D 1.25 sec/record 
Substituting this into (A3) 


Y CF (f) 

y cf^ 


( x 


0.087\ 2 _ 
1.25 / 


0.866 


The measured values of the coherence function, therefore, are 
about 87 percent of the true values. In terms of decibels, then. 
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this represents an error of about 10 log(l/0.866) = 0.6 dB. This 
amount was added to the overall combustor far-f ield power levels 
shown in figure 12 and in the appropriate column of table I, in 
the text. 

Coherence Function Confidence Limits 


Because of the finite sample length of the data, there are 
both variance and bias errors in the computation of the coherence 
functions . 

The confidence limits on the coherence function estimate due 
to these errors may be calculated, to a close degree of approxima- 
tion, from (ref. 38) : 


tanh[w(f) - (n tot - 2 )“H(f) - V= a/2 ] < y^(f) < 

tanh[w(f ) - (n fcot - 2) -1 b(f) + (A6) 

where 

w(f) = tanh -1 [Y CF (f)] 

n. . “ number of statistical degrees of freedom in the estimate = 2N 
tot 


b(f) - [y CF (f) + l.l]/2[y cF (f) + 0.1] 


G w = n/1/ ^ n tot " 2) 

= the 100 a/2 percentage point of the standardised normal 
probability distribution 

Note from the above, that to solve for the confidence limits on 
the coherence function estimate, y^ F (f), the true value, Yq F (£) must 


I 




1 



f 
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be known. The procedure used here to calculate the confidence limits 
was to assume that the estimated value could be substituted for the 
true value in b(£). This is clearly adequate as long as 

tw{f) ± %Z a/2 3(n - 2) » b(f) (A7) 

Furthermore, to produce a conservative result, the estimated 
values of Yqj used in (A6) were those obtained before correcting 
for the time delay bias. (It should be noted that the number of 
statistical degrees of freedom in ref. 38 is exactly one-half the 
number used here. Care must be used, therefore, when comparing 
the equations of ref. 38 with those used here.) 

Example! For 90 percent confidence limits, then 
Z ^/2 = 1.645 (see any table of the standardized normal probability 
distribution) . For the present work with a 120 second taped record, 
and T = 1.25 seconds, then 

N = 120/T = 120 sec/1.25 sec = 96 

n = 2(96) = 192 
tot 

(n tofc - 2)~ 1 = 1/190 = 0.00526 
0 w = n/l/(n fcot - 2) - 0.07267 

For a measured coherence y^(f) of, say, 0.5, we have! 

w(f) - tanh“ 1 (n/0T5) = 0.881 
b(f) - [n/0TI + 1.1]/2[-V0 Ts + 0.1] = 1.119 
Substituting these into (A6), gives: 


tanh[ 0.881 - 0.00526(1.119) - 0.0726(1.645)] < Y cp (f) £ 


or. 


tanh[ 0.881 - 0.00526(1.119) + 0.0726(1.645)] 
0.407 < vL(f) < 0.576 


Or, in decibels, referenced to the estimate 


-0.89 dB < 10 


log 



(f) 

< +0.61 dB 

(f) 


So, for any measured coherence value of 0.5, the associated 
combustor coherence spectrum has a 90 percent confidence limit of 
about -0.9 dB to +0.6 dB just due to the uncertainty in y^(f) . 

For lower values of the measured coherence function, the con- 
fidence limits are wider. For example, at a measured coherence of 
0,1, calculations identical to the above produce 90 percent con- 
fidence limits of -4.1 dB to +2.3 dB. In both cases above, it can 
be verified that the inequality given by (A7) is satisfied. 


Confidence Limits on Combustor Coherence 
Spectra and Power Spectra 

The directivity and power results shown in this study were ob- 
tained by multiplying the measured coherence function by the far- 
field spectrum, on a frequency by frequency basis. The statistical 
errors in the coherence functions were described above. There also 
are, however, statistical errors associated with the far-field spec- 
tral estimate. These errors combine (though not in an additive 
fashion) with the errors in the coherence estimate to produce the 


! 
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final error in the combustor coherence spectra. In reference 26, 
it is shown the variance error, e r> in spectral estimates obtained 
by FFT techniques such as used here is given by: 

M a FF (£)1 - 

where G^Cf) is the far-field spectral estimate. 

For the present case, N => 64 (the far-field spectra were com- 
puted independently, and the number of averages was required to be 
a power of 2 because of the FFT algorithm used by the DSP) , so 

n„ ^ = 128. So, the error is approximately 12.5 percent, or about 
tot 

0.5 dB. This error prevails throughout the spectral range because 
of the constant bandwidth. This is in contrast to an error band of 
at least 1.5 dB associated with the coherence function estimate. In 
reference 30, it is shown that for cases such as this, when one 
error is significantly larger than the other, the larger error dom- 
inates the final coherence spectrum estimate. With the assumption 
that the variance errors in the combustor coherence spectra are dom- 
inated by the variance errors in the coherence function, figure A1 
shows a typical combustor coherence spectrum, with the associated 
90 percent confidence limits sketched in. This figure is for an en- 
gine speed of 43 percent of maximum and a far-field microphone angle 
of 120* . 

With the exception of several frequencies at the 60° far-field 
angle (see appendix C) , the coherence functions between the fluc- 
tuating combustor pressure and the far-field acoustic pressure did 
not vary significantly in either magnitude or shape with respect to 
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far-field angle or engine speed. For this reason, the confidence 
l imi ts on the combustor coherence spectra, as a function of frequen- 
cy, were essentially the same at all far-field angles. The power 
spectra (fig. 53), therefore, which were obtained by a spacial inte- 
gration of the combustor coherence spectra, have 90 percent confi- 
dence limits virtually identical to that of figure Al. The confi- 
dence limits on figure Al, consequently, are representative of all 
spectral results obtained in this investigation via the coherence 
function. 

The bias errors in the far-field spectral estimates are pro- 
portional to the resolution bandwidth, which for this case has an 
effective value of 1 Hz. Because the far-field spectra were all 
smoothly varying (i.e., no tones), these bias errors may be con- 
sidered negligible. 

Confidence Limits on OASPL and OAWPL 

From figure Al, the 90 percent confidence band on the coher- 
ence spectra and the power spectra are quite narrow at frequencies 
near the peak (approximately ±0.5 dB) . This band broadens at fre- 
quencies away from the peak. At 40 Hz the 90 percent confidence 
limits are about -3 dB to +2 dB, and at 240 Hz the 90 percent con- 
fidence limits are about -10 dB to +4 dB. At these two extremes, 
however, the measured spectral levels are considerably below the 
levels at the peak. The large uncertainty band at the extremes of 
the spectra, therefore, do not greatly affect the overall levels. 


which are obtained by anti-logarithmically summing the individual 
spectral levels at each frequency. For example , at the 120° far- 
field microphone angle, for an engine speed of 43 percent of maxi- 
mum, the computed OASFL is 78 dB (fig. 51) , with 90 percent confi- 
dence limits of -1.5 dB to +0.9 dB. Similar results prevail at all 
angles and operating conditions on figure 51. For the OAPWL at 
43 percent operating speed, the computed value is 123.5 dB (fig. 53), 
with 90 percent confidence limits of -1.7 dB to +1.0 dB. Again, 
similar results prevail at the other operating speeds. 

The possibility of random statistical error in the computation 
of the transfer function accounting for the reduced slope cannot be 
estimated since this requires knowledge of the absolute value of the 
measured function (26). Only the relative value is shown in fig- 
ure 21. However, it does not seem likely that random error is re- 
sponsible since all the measured transfer functions with respect to 
the combustor have slopes consistently somewhat less than 2.0. 

Similar comments also apply to the phase results : the confi- 

dence limits on the phase spectra cannot be computed without an ab- 
solute numerical value of the transfer function. However, it is 
true that as the coherence is reduced, the variance of the phase 
estimate increases. This is likely the reason for the apparent 
random phase fluctuation seen at the low and high frequency limits 
of several of the time delay removed phase plots in Chapter 4. 




APPENDIX B 


TABULATION OF COHERENCE FUNCTION CONFIDENCE LIMITS 

This appendix contains a tabulation of the 90 percent confi- 
dence limits for ordinary coherence function estimates between 0.01 
and 0,99. The confidence limits for a specified value of the mea- 
sured coherence are tabulated as a function of the number of inde- 
pendent samples (N) averaged to obtain the estimate. The first set 
of tables is for the linear variable y (f ) . The second set of ta- 
bles is in terms of decibels referenced to the measured value. This 
second set of tables may also be used to determine confidence limits 
on coherence spectra if the variance error in the spectral estimate 
is small compared to the error in the coherence estimate (see appen- 
dix A) . 

Equation (A6) was used to generate the tables for cases where 
the measured coherence is less than 0.75. For tabulated values of 
the measured coherence greater than 0.75, equation (A6) was used with 
b(f) set equal to unity (ref. 38). Blanked locations in the table 
indicate that either N or the measured coherence are too low for 
equation (A6) to be valid. 
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Tabic B-l. 90 Percent Confidence Licice on Coherence Function Measured with Average of N Samples 
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table fl-1. Candy led, 
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Measured 

coherence 

K ■ 

* 16 
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H - 256 
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.650-. 030 

* 6R8-.812 

.712-.799 

+727-. 7BB 

,737-. 760 

+744-. 

775 

,749-. 

770 

+752-.767 

.77 

.599- 

.858 

.663-+BJ7 

.700-. 020 

.723-* 607 

.730-. 797 

.740— .790 

.755-. 

704 

*769-. 

700 

,763-+ 777 

♦ 7B 

.615- 

*B64 

*677-. 345 

. 71 3- • 020 

.735-.B16 

. 749-* 006 

.759-. 799 

.765-. 

794 

+770-+ 

790 

.773-. 731 

.79 

,631- 

.071 

.690-. 052 

* 725-* 836 

*747-* B24 

,760-. H15 

.770- + BUB 

+776-. 

603 

+7B0-. 

799 

+763-.797 

.09 

. 646- 

.877 

.704-. 859 

.730-.B44 

.75B-.033 

+772-. B24 

.700-, P17 

,706-. 

812 

.791-+ 

009 

.793-. 806 

• B l 

. 662- 

.004 

.710-. 067 

*750-. 652 

,770-. 841 

.703-. 833 

+791-+027 

*797-. 

022 

+001-. 

019 

+BO4-.016 

,02 

.679- 

,090 

*732-. B74 

.763- *060 

*702-. 050 

.794-, 042 

+ 602-+83 6 

,000-* 

B31 

.611*. 

028 

+014-.B26 

.63 

,695- 

.097 

*746—* 08 1 

*776-. 666 

.794-.B50 

.BO5-.051 

• 01 3-, 64 5 

.010-. 

041 

.P22-+ 

B3P 

.024-. 036 

.64 

,712- 

.903 

*761-. 0B8 

.7B9-.876 

*006-. B67 

.017-. 860 

.024-. B54 

.829-. 

P50 

.032-. 

HU7 

,935*. 845 

.05 

* 720- 

.909 

.775-. 395 

.001-. 604 

.813-. 075 

. B20-* 86B 

.835-. B63 

-939-. 

860 

.843-. 

957 

+045-.855 

,66 

.745- 

*915 

.7B9-. 903 

.014- * 092 

* 030-* 0B4 

*039-. 877 

+ 04 6- ■ 87 3 

*850-. 

Bfi9 

.053-+ 

066 

.855-, 865 

*67 

, 762- 

*922 

.004-. 910 

, 027- *900 

.041-. 892 

.851-. 886 

+B57-.8R2 

.061-, 

870 

.064-, 

076 

+065-.074 

. 60 

,7 00- 

.928 

. B 18-* 91 7 

. 04O-.9OB 

* 053-+ 901 

.B62-.895 

.06 8-. 09 1 

*871-. 

8BB 

.074-. 

886 

.07&-+0B4 

.69 

.797- 

.934 

.033-. 924 

. 553-, 915 

,066-+ 909 

.873-. 904 

.079-.9OO 

.082-. 

097 

,RP4-+ 

695 

.886-, 094 

,90 

, Q15- 

*940 

.040-. 931 

.866- . 923 

.873-. 917 

.B05-+913 

*090-. 909 

.893-+ 

907 

. 095—. 

905 

,096-. 903 

.91 

. 032- 

*946 

.062-. 93 B 

. 079-. 931 

. 890-.926 

.R96-.921 

.901-.91B 

• 903- • 

916 

• 905*, 

914 

.907-.913 

. 92 

.050- 

.952 

.677-. 945 

.593-. 939 

.902-, 934 

.900-. 930 

.912-.927 

+914—, 

925 

• 9 16—, 

924 

+917-.923 

.93 

.0 60- 

.950 

,392-* 952 

.906-. 946 

.914-. 9«2 

.919-. 939 

.423-. 937 

.925-. 

935 

+926-. 

933 

• 927- .9 32 

. 94 

.087- 

,964 

.900-. 959 

.919-. 954 

.926-, 951 

.931-. 948 

.934-. 446 

.936-* 

944 

*937-. 

943 

. 938-. 942 

.95 

+ 905- 

.970 

.923-. 966 

* 933- * 962 

. 936-. 959 

.942-. 956 

.945-. 955 

-946-+ 

953 

.947-. 

952 

+ 94B-+ 952 

. 96 

.924- 

,976 

.930-. 973 

.946-. 970 

« 951-.967 

*954-. 965 

.956-. 964 

.957-. 

963 

.950-* 

962 

,959-. 961 

.97 

. 943- 

.902 

.953-. 979 

.959-, 977 

,963-. 975 

.965-. 974 

-967-. 973 

.968-. 

972 

.960-+ 

971 

*969-. 971 

. 9B 

,961- 

.900 

. 969— * 906 

*973-. 905 

. 975— *904 

.977-. 903 

. 97B-.9B2 

+976-. 

,961 

. 979-. 

9fl1 

. 979- .961 

.99 

* 9Bl - 

.994 

. 984— ,993 

, 9B6-.992 

*980-. 992 

. 9 B0- + 991 

*989-. 991 

- 909-, 

991 

+969-, 

990 

.990-. 990 
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Table B-2. 90 Percent Confidence Limits on Coherent Output Power Spectrum Computed Vin 

Coherence Function Measured with Average of H Samples (Dccibela) 


Measured 

8 - 16 

N » 32 

N “ 64 

U = LZB 

N - 256 

N - 512 

H “ 1024 

W - 2040 

N “ A 096 

coherence 



















.01 









- 13 , 3 / + 

4 *u 

- 6 . 8 /+ 

3.4 

- 4 * 1 /+ 

2.6 

- 2 . 6 /+ 

1.9 

- 1 . B /+ 

1.4 

.02 









- 6 * B /+ 

3.3 

- 4 . 1 /+ 

2*5 

- 2 . 6 /+ 

1.9 

- 1 . 0 /+ 

1.4 

- 1 . 2 /+ 

1.0 

.03 







- B . 9 /+ 

3.6 

- 5 . 0 /+ 

2.9 

- 3 . 1 /+ 

2.1 

- 2 * 1 /+ 

1.6 

- 1 * 4 /+ 

1.1 

- 1 . 0 /+ 

0.8 

-OB 







- 6 . 9 / + 

3*2 

- 4 . 1 /+ 

2.4 

- 2 , 6 /+ 

1.0 

- 1 . 7 /+ 

1*4 

- 1 * 2 /+ 

1.0 

- U . B /+ 

0,7 

.05 





- 10 * B /+ 

3.7 

- 5 . 7 /+ 

2*9 

- 3 * 5 /+ 

2.2 

- 2 . 3 /+ 

1*6 

- 1 . 5 /+ 

1.2 

- 1 . 0 /+ 

0*9 

- 0 . 7 /+ 

0.6 

.06 





- 8 * 9 / + 

3.4 

- 5 . 0 /+ 

2*7 

- 3 . 1 / + 

2.0 

- 2 * 0 /* 

1.5 

- 1 * 4 /+ 

i.i 

- Q . 9 / + 

0*8 

— 0 . 6 / + 

0,6 

.07 





- 7 . 7 /+ 

3.2 

— 4 , 4 /+ 

2.6 

- 2 . 8 /+ 

1.9 

- 1 * 3 /+ 

1.4 

- 1 . 2 / + 

1.0 

- 0 . 8 /+ 

0*7 

- 0 . 6 /+ 

0.5 

.06 





- 6 , B /+ 

3.0 

- 4 * 0 /+ 

2.3 

- 2 , 5 /+ 

1*7 

— 1 * 7 /+ 

1.3 

- 1 * 1 /+ 

0.9 

- O . 0 /+ 

0.7 

- 0 . 5 /+ 

0.5 

♦ 09 



- 12 . 4 / + 

3*6 

- 6 . 1 /+ 

2 . B 

- 3 . 7 /+ 

2.2 

- 2 . 3 /+ 

1.6 

- 1 . 6 /+ 

1.2 

- 1 . 1 /+ 

0*9 

- 0 . 7 /+ 

0*6 

- Q . 5 /+ 

0*5 

.10 



- 10 . 8 /+ 

3.4 

- 5 . 6 /+ 

2*7 

- 3 * 4 /+ 

2.0 

- 2 . 2 /+ 

1*5 

- 1 * 4 /+ 

1.1 

- 1 . 0 /+ 

O.B 

- 0 . 7 /+ 

0.6 

- 0 . 5 /+ 

0.4 

*11 



- 9 * 7 /+ 

3.3 

- 5 . 2 /f 

2.5 

- 3 , 2 /+ 

1*9 

- 2 . 0 /+ 

1,5 

- 1 . 4 /+ 

1.1 

- 0 . 9 /+ 

0*8 

- 0 . 6 /+ 

0.6 

—0. U/+ 

0,4 

.12 



- B . 8 /+ 

3.1 

- 4 . B /* 

2.4 

- 3 * 0 /+ 

1.8 

- 1 . 9 /+ 

1.4 

- 1 . 3 /+ 

1.0 

- Q . 9 /+ 

0.7 

- 0 . 6 /+ 

0.5 

- 0 . 4 /+ 

0*4 

.13 



- 8 . 1 /+ 

3.0 

- 4 * 5 /+ 

2.3 

- 2 . B /+ 

1*0 

- 1 . 8 /+ 

1.3 

- 1 * 2 /+ 

1.0 

- 0 . 8 /+ 

0.7 

- 0 , 6 /+ 

0.5 

- 0 . 4 /+ 

D .4 

. m 



- 7 . 5 /+ 

2.9 

- 4 . 2 /+ 

2.2 

- 2 * 6 /+ 

1.7 

- 1 . 7 /+ 

1.3 

- 1 . 2 /+ 

0.9 

- 0 . 8 /+ 

0.7 

- 0 . 5 /+ 

0.5 

— 0 , 4 /+ 

0.3 

.15 

^ 16 . 6 /+ 

3.5 

- 7 . 0 /+ 

2.8 

- 4 . 0 /+ 

2.1 

- 2 . 5 /+ 

1.6 

-1 . 6 /+ 

1.2 

- 1 . 1 /+ 

0.9 

- 0 - 7 /+ 

0.6 

- 0 . S /+ 

0.5 

- 0 . 4 /+ 

0*3 

.16 

- 14 . 0 /+ 

3.3 

- 6 . 5 /+ 

2,7 

- 3 . 8 /+ 

2 . 1 

- 2 . 4 /+ 

1.6 

- 1 . 6 / + 

1.2 

- 1 . 0 /+ 

0.9 

- 0 . 7 /+ 

0*6 

- 0 . 5 /+ 

0*4 

- 0 * 3 /+ 

0.3 

.17 

- 13 . 4 /+ 

3.2 

- 6 . 2 /+ 

2.6 

- 3 . 6 /+ 

2.0 

- 2 . 3 /+ 

1,5 

- 1 . 5 /+ 

1.1 

- 1 . 0 /+ 

0.0 

- 0 . 7 /+ 

0*6 

- 0 . 5 /+ 

0*4 

- 0 . 3 /+ 

0.3 

.16 

- 12 . 2 /* 

3*1 

- 5 . B /+ 

2.5 

- 3 . 4 /* 

1*9 

- 2 . 2 /+ 

1.4 

- 1 . 4 /+ 

1.1 

- 1 . 0 /+ 

O.B 

- 0 . 7 /+ 

0,6 

- 0 * S /+ 

0.4 

- 0 . 3 /+ 

0.3 

.19 

- 11 , 3 /+ 

3*0 

- 5 . 5 /t 

2.4 

- 3 . 3 / + 

1 .B 

- 2 * 1 / + 

1.4 

- 1 . 4 / + 

1,0 

- 0 . 9 /+ 

0.0 

- 0 . 6 /+ 

0.5 

- 0 * 4 /+ 

0.4 

- 0 , 3 /+ 

0.3 

.20 

- 10 . 5 /+ 

2.9 

- 5 . 2 / + 

2.3 

- 3 * 1 /+ 

1.8 

- 2 . 0 /+ 

1.3 

- 1 * 3 /+ 

1.0 

- 0 . 9 /+ 

0.7 

- 0 * 6 /+ 

0.5 

- 0 . 4 / + 

0.4 

-0* v* 

0.3 

*21 

- 9 . 9 /+ 

2*8 

- 5 . 0 /+ 

2.2 

- 3 . 9 /+ 

1.7 

- 1 . 9 /+ 

1*3 

- 1 . 3 / + 

1*0 

- 0 . 9 /+ 

0.7 

- 0 . 6 /+ 

0.5 

- 0 . 4 /+ 

0.4 

- 0 . 3 /+ 

0.3 

.22 

- 9 . 3 /+ 

2*7 

— 4 . 6 /+ 

2.2 

- 2 , 9 /+ 

1,7 

- 1 * 8 /+ 

1.3 

- 1 . 2 /+ 

0*9 

- 0 * 0 /+ 

0.7 

- 0 . 6 /+ 

0,5 

- 0 * 4 /+ 

0.4 

- 0 * 3 /+ 

0*3 

*23 

- 0 . 8 /* 

2.7 

- 4 . 6 /+ 

2.1 

- 2 * 7 /+ 

1.6 

- 1 . B /+ 

1.2 

- 1 . 2 / + 

0.9 

- 0 . B /+ 

0.7 

- 0 * 5 /+ 

0.5 

- 0 . B /+ 

0.3 

- 0 * 3 / + 

0.2 

.24 

- 8 . 3 /* 

2,6 

- 4 * 4 / + 

2.0 

- 2 . 6 /+ 

1.6 

- 1 * 7 /+ 

1.2 

- 1 . 1 / + 

0.9 

- 0 . B /+ 

0.6 

- 0 . 5 /+ 

0.5 

- 0 . 4 / + 

0,3 

-0*V+ 

0.2 

.25 

- 7 . 9 /+ 

2.5 

- 4 . 2 /+ 

2.0 

- 2 . 5 /+ 

1.5 

- 1 . 6 /+ 

1.1 

-i. i/+ 

0 .B 

- 0 . 7 /+ 

0.6 

- 0 . 5 /+ 

0.4 

- 0 . 4 /+ 

0,3 

- 0 . 2 /+ 

0*2 

.26 

- 7 . 5 /+ 

2,4 

- 4 . 0 /+ 

1.9 

- 2 * 4 / + 

1 . 5 

- 1 * 6 /+ 

1*1 

- 1 * 0 /+ 

0.8 

- 0 - 7 /+ 

0.6 

- 0 . 5 /+ 

0.4 

- 0 . 3 /+ 

0*3 

- 0 . 2 /+ 

0.2 

.27 

- 7 . 2 /+ 

2.4 

- 3 . 9 /+ 

1.9 

- 2 * 4 /+ 

1*4 

- 1 . 5 /+ 

1.1 

- 1 . 0 /+ 

0,8 

- 0 . 7 /+ 

0.6 

- 0 . 5 /+ 

0*4 

- 0 * 3 /+ 

0.3 

- Q . 2 /+ 

0*2 

* 2 b 

- 6 . B /+ 

2.3 

- 3 . 7 / + 

1.8 

- 2 - 3 /+ 

1.4 

- 1 . 5 /+ 

1.0 

- 1 . D /+ 

0.3 

- 0 . 7 /+ 

0.6 

- 0 * 5 /+ 

0.4 

- 0 . 3 /+ 

0.3 

- 0 * 2 /+ 

0.2 

.29 

— 6 * 6 /* 

2*2 

- 3 * 6 /+ 

1.1 

- 2 , 2 /* 

1*3 

- 1 * 4 /+ 

1*0 

- 0 . 9 /+ 

0*7 

- 0 . 6 /+ 

0.5 

- 0 . 4 /+ 

0.4 

- 0 . 3 /+ 

0*3 

- 0 . 2 /+ 

0*2 

.30 

- 6 . 3 /+ 

2*2 

— 3 . 5 / + 

1.7 

- 2 . 1 /+ 

1*3 

- 1 * 4 /+ 

1.0 

- 0 . 9 /+ 

0.7 

- 0 . 6 /+ 

0.5 

-Q.4/+ 

0.4 

- 0 * 3 /+ 

0 * 3 

- C . 2 /+ 

0*2 

.31 

- 6 . 0 /* 

2 * 1 

- 3 . 3 /+ 

1.7 

- 2 * 1 /+ 

1*3 

- 1 . 3 /+ 

0.9 

- 0 . 9 / + 

0,7 

- 0 * 6 /+ 

0*5 

- Q . 4 /+ 

0.4 

- 0 . 3 /+ 

0.3 

- 0 * 2 /+ 

0.2 

.32 

- 5 * 8 /+ 

2.0 

- 3 * 2 /+ 

1.6 

- 2 * 0 /+ 

1.2 

- 1 , 3 /* 

0,9 

- 0 . 9 /+ 

0.7 

- 0 . 6 /+ 

0*5 

- 0 * 4 /+ 

0*4 

- 0 . 3 /+ 

0.3 

- 0 * 2 /+ 

0*2 

*33 

- 5 . 6 /+ 

2.0 

- 3 . 1 /+ 

1.6 

- 1 . 9 /+ 

1.2 

- 1 * 2 /+ 

0,9 

- 0 . B /+ 

0.7 

- 0 . 6 /+ 

0*5 

- 0 * 4 /+ 

0.3 

- 0 . 3 /+ 

0*2 

- 0 . 2 /+ 

0.2 

. 34 

- 5 . 4 /+ 

7*9 

- 3 . 0 /+ 

1.5 

- 1 . 9 /+ 

1.2 

- 1 * 2 /+ 

0.9 

- 0 . B /+ 

0.6 

- 0 . 6 /+ 

0*5 

- 0 * 4 /+ 

0.3 

- 0 - 3 /+ 

0*2 

- 0 * 2 /+ 

0*2 

.35 

- 5 . 2 /* 

1.9 

- 2 . 9 /* 

1.5 

- 1 - B /+ 

1,1 

- 1 . 2 /+ 

0*8 

- 0 . B /+ 

0,6 

- 0 . 5 /+ 

0.5 

- 0 . 4 /+ 

0*3 

- 0 . 3 /+ 

0*2 

- 0 . 2 /+ 

0*2 

.36 

- 5 . 0 /+ 

1.0 

- 2 - B /+ 

1.4 

- 1 . 7 /+ 

1.1 

- 1 * 1 /+ 

0.0 

- 0 . 8 /+ 

0.6 

- 0 . 5 /+ 

0*4 

- 0 . 4 /+ 

0.3 

- 0 . 2 /+ 

0*2 

- 0 . 2 /+ 

0.2 

.37 

- 4 . B /+ 

1.8 

- 2 * 7 /+ 

1.4 

- 1 . 7 /+ 

1,1 

- 1 * v + 

0.8 

- 0 . 7 /+ 

0*6 

- 0 . 5 /+ 

0*4 

- 0 . 3 /+ 

0.3 

- 0 * 2 /+ 

0*2 

- 0 , 2 /+ 

0.2 

.30 

- 4 . 6 /+ 

1.7 

- 2 . 6 /+ 

1.4 

- 1 . 6 /+ 

1.0 

- 1 . 1 /+ 

0 .B 

- 0 . 7 /+ 

0*6 

- 0 . 5 /+ 

0,4 

- 0 . 3 /+ 

0.3 

- 0 . 2 /* 

0*2 

- 0 * 2 /+ 

0.2 

,39 

- 4 * 5 /+ 

1.7 

- 2 . 6 /+ 

1,3 

- 1 , 6 / + 

1*0 

- 1 . 0 /+ 

0.3 

- 0 , 7 / + 

0,6 

- 0 . 5 /+ 

0.4 

- 0 * 3 /+ 

0.3 

- 0 . 2 /+ 

0.2 

- 0 - 2 /+ 

0.1 

.40 

- 4 . 3 /+ 

1-6 

- 2 * 5 /+ 

1.3 

- 1 * 5 /+ 

1*0 

- 1 * 0 /+ 

0*7 

- 0 * 7 /+ 

0.5 

- 0 . 5 /+ 

0*4 

- 0 . 3 /+ 

0.3 

- 0 . 2 /+ 

0.2 

- D . 2 /+ 

C ,1 

.41 

- 4 . 2 /+ 

1.6 

— 2 * 4 /+ 

1.3 

- 1 * 5 /+ 

1.0 

- 1 * 0 /+ 

0.7 

- 0 . 7 /+ 

0.5 

- 0 * 4 /+ 

0*4 

- 0 . 3 /+ 

a. 3 

- 0 * 2 /+ 

0.2 

- 0 . 1 /+ 

0.1 

.42 

- 4 . 0 /+ 

1.6 

- 2 . 3 / + 

1.2 

- 1 * 4 /+ 

0*9 

- 0 * 9 /+ 

0.7 

- 0 * 6 /+ 

0*5 

- 0 . 4 /+ 

0*4 

- 0 . 3 /+ 

a, 3 

- 0 . 2 /+ 

0,2 

- 0 . 1 /+ 

0 . 1 

.43 

- 3 . 9 / + 

1.5 

- 2 . 2 /+ 

1*2 

- 1 . 4 /+ 

0.9 

- 0 . 9 / + 

0.7 

- 0 , 6 /+ 

0*5 

- 0 . 4 /+ 

0 - 4 

- 0 * 3 /+ 

0 . 3 

- 0 * 2 /+ 

0*2 

- 0 * 1 /+ 

0.1 

.44 

- 3 . 8 /+ 

1.5 

- 2 . 1 /+ 

1*2 

- 1 . Q /+ 

0.9 

- 0 , 9 /+ 

0*7 

- 0 . 6 /+ 

0.5 

- 0 . 4 /+ 

0,4 

- 0 * 3 /+ 

0.3 

- 0 . 2 /+ 

0,2 

- 0 * 1 /+ 

0*1 

.45 

-3.6/+ 

1*4 

- 2 . 1 /+ 

1.1 

- 1 . 3 /+ 

0*9 

- 0 . 9 /+ 

0*6 

- 0 . 6 /+ 

0.5 

- 0 * 4 /+ 

0.3 

- 0 * 3 /+ 

0.2 

- 0 . 2 /+ 

0.2 

- 0 . 1 /♦ 

0.1 

.46 

- 3 . 5 /* 

1.4 

- 2 . 0 /+ 

1.1 

- 1 * 3 /+ 

0,8 

- 0 . 8 /+ 

0.6 

- 0 . 6 /+ 

0.5 

- 0 . 4 /+ 

0.3 

- 0 . 3 /+ 

0.2 

— 0 + 2 /+ 

0.2 

-0.1/+ 

0.1 

.47 

- 3 . 4 /+ 

1,4 

- 2 . 0 /+ 

1.1 

- 1 . 2 /* 

0.0 

- 0 . B /+ 

0.6 

- 0 * 5 /+ 

0.4 

- 0 . 4 /+ 

0.3 

- 0 . 3 /+ 

0.2 

- D . 2 /+ 

0.2 

- 0 . 1 /* 

0.1 

.48 

- 3 . 3 /* 

1*3 

- 1 . 9 /+ 

1*0 

- 1 . 2/4 

0.0 

- 0 * 8 /+ 

0.6 

- 0 , 5 /+ 

0.4 

- 0 . 4 /+ 

0.3 

- 0 . 3 /+ 

0.2 

- 0 * 2 /+ 

0.2 

-0.V+ 

0.1 

.49 

- 3 . 2 /+ 

1-3 

- 1 . 0 /+ 

1*0 

- 1 * 2 /+ 

0.3 

- 0 . 0 /+ 

0.6 

- 0 * 5 /+ 

0.4 

- D - 4 /+ 

0.3 

- 0 * 2 /+ 

0.2 

- 0 . 2 /+ 

0.2 

- 0 . 1 /+ 

a.i 

.50 

- 3 . 1 /+ 

1.2 

- 1 . 0 /+ 

1*0 

- 1 . 1 /+ 

0*7 

- Q . 7 /+ 

0.6 

- 0 * 5 /+ 

0.4 

- 0 * 3 /+ 

0*3 

- 0 . 2 /+ 

0*2 

- 0 , 2 /+ 

0.2 

- 0 . 1 /+ 

0*1 
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table B-2. Concluded, 


Measured 

coherence 

H - 16 

N * 32 

H - 64 N - 128 

N - 256 

U - 512 

N = 1024 

N * 2040 

N - 4096 

*51 

- 3 , 0 /+ 

1.2 

- 1 . 7 /+ 0.9 

- 1 . 1 /+ 

0.7 - 0 . 7 /+ 

0.5 

- 0 . 5 /+ 

0.4 

- 0 . 3 /+ 

0.3 

- 0 * 2 /+ 

0.2 

- 0 , 2 /+ 

0-1 

- 0 . 1 /+ 

0 . 1 

*52 

- 2 . 9 /+ 

1.2 

- 1 * 7 /+ 0.9 

- 1 * 1 /+ 

0*7 - 0 . 7 /t 

0.5 

- 0 . 6 /+ 

0*4 

- 0 . 3 /+ 

0*3 

- 0 . 2 /+ 

0*2 

- 0 . 2 /+ 

0*1 

- 0 . 1 /+ 

0.1 

.53 

- 2 . e /* 
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APPENDIX C 




DIRECTIVITY PATTERN FOR SEVERAL FREQUENCIES 

In chapter 5 it was shewn that the OASPL combustion directivit- 
ies as obtained from the coherence functions agreed reasonably well 
with those of most previous studies. The directivity patterns were 
relatively uniform, with a slight peak near 120° (relative to the 
engine inlet axis) . The same was not true, however, for certain 
frequencies. The coherence function at 43 percent engine speed be- 
tween combustor pressure and the far-field acoustic pressure at the 
60° microphone angle was shown in figure 36 and was seen to have 
several pronounced dips at frequencies of 40, 80, 120, and 160 Hz. 

The corresponding coherence spectrum was shown in figure 50(a) and, 
because of the relatively flat far-field spectrum at this angle, also 
contained dips at these frequencies. These nips gradually disap- 
peared at angles away from 60°. The result, however, is an unusual 
directivity pattern at these frequencies, and this is shown in fig- 
ure Cl. In addition to the peak near 120° , there is now a 7 to 
10 dB drop at 60°. Similar results were found for these same fre- 
quencies at all engine operating speeds through 60 percent. This 
dip does not show up in the OASPL directivity, of course, because 
of the logarithmic addition associated with obtaining an OASPL. 

There are several possible reasons for this result. The first 
of these simply involves the presence of a noise source at these 
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frequencies which does not correlate with the combustor pressures. 
Such an extraneous noise source, however, would have to be highly 
directional (l.e., radiating strongly only at the 60* angle) and 
this seems quite inconsistent with the low frequencies involved. 

It is also possible that the combustion noise is radiating from 
both aft and front ends of the engine, arriving at the 60* micro- 
phone out of phase, resulting in a partial cancellation. This satu- 
ration may arise from the time delay difference between the forward 
traveling acoustic waves (the speed of which is diminished by the 
gas flow speed within the engine) and the rearward traveling acoustic 
waves (the speed of which is augmented by the gas flow speed within 
the engine) . Additionally, the plane of the bellmouth inlet is at 
a different distance from the 60° microphone than is the exit plane 
of the primary nozzle. The combined difference in propagation time 
may be sufficient to introduce a phase shift at the 60° microphone 
which results in cancellation. Such a situation, however, would ob- 
viously have to prevail only at the frequencies in question. 

Finally, we are left with the possibility of some unusual phas- 
ing relationship at 60® between casing radiation (i.e., the engine 
as a unit vibrating) and direct acoustic radiation from the aft end. 
The highest frequency in question (200 Hz) has a corresponding wave- 
length of about 1.7 meters. This is certainly not large in compari- 
son with the length of the engine, which is approximately 3.75 me- 
ters long. Phasing between casing radiation and aft-end radiation 
to produce a partial cancellation of 60° is, therefore, quite pos- 
sible. Again, however, this phenomenon would have to be restricted 
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only to those frequencies in question. 

The final remark to be made in this appendix, therefore, is 
simply that no definite conclusions can be reached regarding the dip 
in the directivity pattern which occur only at very specific fre- 
quencies, independent of engine operating speed. Precisely the same 
results were observed when the tests were rerun about 10 days later. * 

The possibility, therefore, of an isolated anomalous test arena or 
environmental condition may be ruled out. 




APPENDIX D 


FOURIER ANALYSIS OF RANDOM SIGNALS 

The purpose of this appendix is to provide a brief review of 
the various frequency and time domain representations of random sig- 
nals. It is intended primarily as background information for the 
reader not familiar with the techniques of random data analysis. 
Rigorous derivations and detailed interrelationships of the various 
functions used in the main body of this thesis, including digital 
computation techniques, may be found in numerous textbooks on the 
subject. Reference (26) is the primary source for most of the mate- 
rial presented in this appendix. The Fourier transform operation 
referred to here is understood to be the discrete, finite range 
Fourier transform described, for example, in ref. (26). 

Auto-Spectral Density 

Given a single realization of a random signal, say x(t) , 
then its Fourier transform S (jw) will, in general, have a real 
and imaginary part: 

S (jto) = -$[x(t)] = a (to) + jb(to) (j = ~/~l) 

where to is the circular frequency. 

The auto-spectral density of the signal G (to) may be approxi- 


mated by: 


G (w) “ 55 a2 (“) + b2 (w) 

•tUv 3C 

where the * signifies complex conjugate. 

The auto-spectral density is a real-valued function, and is the 
quantity measured by conventional analog filtering procedures, as 
the filter bandwidth becomes very small. 

The basic method for implementing this calculation digitally is 
as follows: 

1. Record a sample of the signal 

2. Do a Fourier transform of the data sample 

3. Multiply the Fourier transform by its complex conjugate 

4. Add this result to the previous sum of spectral densities 

5. Repeat the sequence above N times 

6. Divide the final result by N to obtain an ensemble averaged 

spectral density estimate 

Note that because the signal is random, an estimate of the true spec- 
tral density is obtained by averaging the results of N ensembles. 
The spectral densities themselves must be averaged and not the 
Fourier transforms, since the real and imaginary parts would average 
to zero for a random signal, 

Cross-Spectral Density 

The cross-spectral density G^Cjui) between two signals, say 
x(t) and y(t), indicates similarities between the two signals in 
the frequency domain and contains relative phase information. It is 
obtained from: 
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G^Cjoi) = S x (ju)S*(ja)) 

In general, this will be a complex quantity with real and imag- 
inary parts, sometimes referred to as the coincident and quadrature 
spectral densities, respectively. That is 

G^Cjio) = Co(w) + j Quad (to) 

The magnitude of the cross-spectral density is real valued and 
for random functions is given by the polar combination of the aver- 
aged Co and Quad spectra: 


iG^CjoOl 


[Co 2 (to) + Quad 2 (to) ] 


1/2 


The phase of the cross -spectral density is also real-valued and is 
given by 


6 Oo) “ 0„(m) - 6 Co) =■ Arctan 
xy x y 


Quad(m) 


i 


Co (id) 

The phase of the cross-spectral density is the phase difference be- 
tween the signals x(t) and y(t). 


Transfer Function 


The transfer function H (jw) between two signals x(t) and 

xy 

y(t) provides frequency domain information about the nature of the 
transfer of a signal through a system. It Is a -complex quantity hav- 
ing real and imaginary parts and is given by the cross spectral den- 
sity between the two signals normalized with respect to the auto- 

spectral density of the signal considered as input: 

G_Cju) 

V (:|W) ' G^Cw) 
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Its magnitude, |h (joi) j , is real valued and is given by polar 
xy 

combination of the real and imaginary parts. Its phase is also real 
valued and given by the inverse tangent of the ratio of the imaginary 
to real parts, and is clearly the same as the phase of the cross- 
spectral density. 


Ordinary Coherence Function 


The ordinary coherence function y (ai) between two random sig- 

xy 

nals, say x(t) and y(t), is a numerical measure of the similari- 
ties between the two signals at each frequency. It is given by the 
square of the ensemble averaged cross-spectral density normalized 
with respect to the product of the individual averaged auto-spectral 
densities . 

2 l°„W | 2 

It is a real valued function which must take on a value between 
zero and one, with high (low) coherence at a particular frequency 
meaning high (low) correlation at that frequency. As described in 
chapter 5 it may, under the appropriate circumstances, be inter- 
preted as the fractional portion of the mean-square value at y(t) 
which is due to the input x(t) . It should be pointed out that when 
computing coherence functions from spectral density estimates, as 
defined in the equation above, the ensemble averaging must be per- 
formed on the individual spectral densities, and not on the computed 
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coherence function from each sample. The latter approach will al- 
ways yield a coherence of unit (see ref. (26)). 

Cross-Correlation Function 

The time domain definition of the cross-correlation function, 
for random signals is: 

R (t) = lin x(t)y(t + x)dt 

The cross-correlation function and the cross-spectral density 

are Fourier transform pairs and hence R (t) may also be computed 

xy 

from: 

The DSP used in this investigation utilized the Fourier 
transform relationship above to compute the cross-correlation 
functions . 


b 


I 
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